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Power Corrections to Event Shapes with Mass-Dependent Operators 
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We introduce an operator depending on the "transverse velocity" r that describes the effect of 
hadron masses on the leading 1/Q power correction to event-shape observables. Here, Q is the scale 
of the hard collision. This work builds on earlier studies of mass effects by Salam and Wicke [lj and of 
operators by Lee and Sterman Q| . Despite the fact that different event shapes have different hadron 
mass dependence, we provide a simple method to identify universality classes of event shapes whose 
power corrections depend on a common nonperturbative parameter. We also develop an operator 
basis to show that at a fixed value of Q, the power corrections for many classic observables can be 
determined by two independent nonperturbative matrix elements at the f 0% level. We compute the 
anomalous dimension of the transverse velocity operator, which is multiplicative in r and causes the 
power correction to exhibit non-trivial dependence on Q. The existence of universality classes and 
the relevance of anomalous dimensions are reproduced by the hadronization models in Pythia 8 and 
Herwig++, though the two programs differ in the values of their low-energy matrix elements. 



I. INTRODUCTION 

Event shapes have played a key role in establishing the 
structure of quantum chromodynamics (QCD). Indeed, 
the jet-like behavior of QCD collisions at high energies 
was established in 1975 by measuring the event shape 
sphericity in e + e~ collisions Q. Event shapes are used 
for precision determinations of the strong coupling con- 
stant a s using e + e~ data Q, and the same event shapes 
are used to tune Monte Carlo hadronization models to 
help make particle-level predictions for the Large Hadron 
Collider (LHC) (see e.g. Ref. @). Recently there has also 
been a resurgence of interest in event shapes because they 
are closely related to jet shapes, which are a sensitive 
probe of jet substructure @, . 

Event shapes require both perturbativc and nonpertur- 
bative contributions from QCD. Therefore power correc- 
tions are an important theoretical ingredient for making 
event shape predictions. In the kinematic tail region, the 
leading nonperturbative effect is simply to shift perturba- 
tive event shape distributions by an amount suppressed 
by 1/Q [8,]. Such 1/Q n effects are known as power cor- 
rections, where Q is the scale of the hard collision. Some 
power corrections are known to exhibit universality, in 
the sense that the same nonperturbative shift parameter 
describes more than one event shape. Power corrections 
also encode the effect of hadron masses on event shape 
distributions since different methods for treating the 
energy /momentum of soft hadrons with mass mu can 
change event shapes by 0(m,H /Q)- 

In this paper, we revisit the theoretical underpinnings 
of power corrections for e + e~ event shapes, with a partic- 
ular focus on universality and hadron masses. Building 
on the work of Salam and Wicke [l[ , we show that hadron 
mass effects break traditional power correction universal- 
ity, but we are still able to define universality classes of 
event shapes with a common power correction. Building 
on the work of Lee and Sterman Q , we provide an opera- 
tor definition of the leading power correction that demon- 



strates that its hadron mass dependence is described by 
a "transverse velocity" distribution. Transverse velocity 
r is defined as 
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where the transverse momentum p± is measured with re- 
spect to the thrust axis. The transverse velocity operator 
has nontrivial renormalization group evolution which is 
multiplicative in r, such that the leading power correc- 
tion includes a term behaving as (a s h\Q)/Q. 

To put our work in context, it is worthwhile to review 
some of the key literature on event shapes and power 
corrections. Examples of classic e + e~ dijet event shapes 
arc thrust j^L the C-paramete r fiol llil ] , hemisphere jet 
masses fl2 - fl4l ] , jet broadening [15| , and angularities [16| , 
all of which were measured at LEP. The theoretical un- 
derstanding of event shapes has undergone substantial 
progress in recent years, especially for perturbative cal- 
culations. Fixed-order corrections to event-shape distri- 
butions have been calculated up to 0(cv 3 ) precision 
[23j , and large singular logs have been resummed to N 3 LL 
accuracy for thrust [24| and heavy jet mass [25|, and to 
NLL for jet broadening p6l - [28| . For certain event shapes 
e, it can be demonstrated that the leading power correc- 
tion parameter f2f appears both in the mean value of the 
event shapes as well as in the tail region of the full event 
shape distribution This fact was discussed using a 
factorization theorem for thrust in Ref. [2^] , and used to 
simultaneously extract a s (mz) and f2[ from thrust data 
at N 3 LL + e>(a 3 ) in Refs. HI 

To study nonperturbative power corrections, there 
have been two broad strategies. The first strategy is 
to build analytic models to describe the nonperturbative 
physics. For example, in the dispersive approach (3~l| - 
|33[, one introduces an infrared cutoff fj,j below which 
the strong coupling constant is replaced by an effective 
coupling ev ff. Perturbative infrared effects coming from 
scales below the cutoff are subtracted, and nonperturba- 
tive effects are parametrized in terms of an average value 
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for the effective coupling ao times /^/Q The original dis- 
persive approach relied on the single-gluon approxima- 
tion. In Refs. [HI, H3] the Milan factor was introduced 
as a refinement to handle multi-gluon diagrams. An- 
other renormalon-inspired a naly tic model is provided in 
the dressed gluon approach [4ll - l43l |. 

A key prediction that emerged from these analytic 
models (first seen in the dispersive approach) was univer- 
sality of power corrections [3l|, [35| (see also [39[ ) . Uni- 
versality posits that the leading power correction flf to 
an event shape e can be separated into two pieces, a cal- 
culable coefficient c e which depends on the event shape 
in question times a universal nonperturbative parameter 
A common to all event shape observables: 

n\ -> c e A. (2) 

Because of universality, one can extract the power cor- 
rection A from one event shape measurement and apply 
it to another one, as studied in Refs. [2I I3I liol. lima] . 

The second strategy to understand power corrections 
is based on QCD factorization, where there is a formal 
separation of pcrturbativc and nonperturbative contribu- 
tions. The shape function was introduced in Refs. [47l.l49j 
to describe nonperturbative corrections, and it accounts 
for a whole range of power corrections of the order 
(AQCD/eQ)™. In the tail region where Aqcd < Qe < Q, 
the shape function can be expanded in terms of deriva- 
tives of the Dirac delta function, and this translates into 
an operator expansion for the event-shape distribution. 
This shape function can be derived in the Collins- Soper- 
Stcrman (CSS) approach to factorization [l|| H3, l49l - 
[5l| . The shape function also emerges naturally from fac- 
torizatio n prope rties of Soft-Collinear Effective Theory 
(SCET) p52T - l56| . Here methods exist to systemati cally 
improve the description of the shape function [57], l58| . 
and to sum large logs present in the subtractions needed 
to define the power corrections in a scheme that is free 
from the leading renormalon ambiguity [59l |60| . 

The key advantage of the factorization-based strategy 
is that one can express nonperturbative parameters in 
terms of matrix elements involving QCD fields. For ex- 
ample, the leading power correction fix can be expressed 
in terms of the energy flow operator 0, loH - 1631 ]. Using 
this fact, Lee and Sterman Q rigorously derived univer- 
sality for the leading power correction without relying 
on an analytic model. This proof of power correction 
universality solely uses QCD first principles. 

The drawback of both of the two above strategies is 
that they start by making the assumption that all final 



The dispersive approach was motivated by the study of renor- 
malons (see Ref. |34| for a review). Renormalons refer to an 
ambiguity in a resummed perturbative series of order (A/Q) p , 
and this ambiguity is related to the nonperturbative power cor- 
rection. Rcnormalon-bascd models were originally applied to the 
mean values of the event-shape distributions 31, 35- 37] (sec also 
f3Sj ~) and later generalized to the dijet limit of distributions as 

well Hm. 



state hadrons can be considered to be massless. An im- 
portant exception to this is Ref. [13, H3, where mass 
effects were studied with renormalons using massive glu- 
ons. This is certainly true at the parton level, but the 
actual hadrons measured in an event shape are of course 
massive. One might erroneously argue that at very high 
energies the hadron masses can be safely neglected, but 
because nonperturbative effects are caused by soft par- 
ticles whose momenta are of order Aqcd ~ "iff, their 
masses contribute at order mn/Q, which is the same or- 
der as the leading power correction. One could try to 
study the effect of hadron masses by using Monte Carlo 
hadronization models, but this approach is not fully sat- 
isfactory since there is no clean separation between per- 
turbative parton shower evolution and nonperturbative 
hadronization effects, and there is no guarantee that one 
can systematically improve the accuracy of the treatment 
of hadronization effects. 

The first serious study of hadron mass effects on power 
corrections was performed by Salam and Wicke [l| , where 
they argued that universality does not hold for event 
shapes such as thrust, jet masses, or the C-parameter 
with their traditional definitions. However, they showed 
that universality can be restored if one makes measure- 
ments in the E-scheme, where one performs the following 
substitution in the event shape definition: 

E 

Pi 7^7 P* ■ (3) 
\Pi\ 

Ref. [T| also argues that hadron mass effects, in addition 
to breaking universality, generate a power correction of 
the form (\nQ) A /Q, with A ~ 1.5. 

In this paper we will devise a rigorous operator-based 
method of treating hadron mass effects in the leading 
event shape power correction by defining a transverse ve- 
locity operator. We generalize fif to a function Oi(r, fi) 
which accounts for both the transverse velocity depen- 
dence through r as well as renormalization group evolu- 
tion through /1. Using this operator we derive univer- 
sality classes for event shapes in the presence of hadron 
masses, following the treatment of Ref. 0, and show that 
each event shape belongs to a unique class. From study- 
ing the r- and /x-dependence, we will largely confirm the 
results of Salam and Wicke. 

The remainder of this paper is organized as follows. In 
Sec. HIl we set the notation for event shapes and power 
corrections. In Sec. IIII[ we introduce the transverse ve- 
locity operator whose matrix elements yield the leading 
power correction. We explore the consequences of univer- 
sality classes in Sec. lIVi and use a complete operator basis 
to derive approximate universality relations. We consider 
the effect of renormalization group evolution in Sec. fVl 
where we derive the anomalous dimension for fii(r, fi) 
and compare to the hadronization models of Pythia 8 
and Hcrwig++. We conclude in Sec. IVII with a discus- 
sion of the implications and extensions of our results, in 
particular for the LHC. 
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Thrust 2 : 


r = 






Qp T 

T= Q = 




2-jettiness Hi]: 


T~2 = 






T2 = T 2 = 




Angularities 


T(a) = 




- Icos^l) 1 -- , 


f(a) = r (o) = 


1 P 


C-paramcter [l^, EH 3 : 


c = 


p i,j 




c = 


3 

Q ^ cosh(r/ l ) ' 


Jet Masses (TMl: 


p± = 






p± = 





TABLE I. Examples of event shapes e with their original definitions, as well as formulae e that are valid for e -C e max . See 
the text below Eq. (0 for a further description of the notation. 



II. POWER CORRECTIONS FOR DUET 
EVENT SHAPES 

We begin by reviewing the notation for kinematics and 
event shapes in Sec. Ill Al various hadron mass schemes in 
Sec. Ill Bl and the basics of how power corrections impact 
event shapes in Sec. Ill Cl Readers familiar with these top- 
ics can skip to Sec. lIIII where we introduce the transverse 
velocity operator. 



A. Event Shapes with Transverse Velocity 

A dimensionless event shape e is an observable defined 
on final-state particles which can be used to describe the 
jet-like structure of an event. To describe particle mo- 
menta, we use rapidity y, pseudo-rapidity 77, transverse 
momenta p± = \p±\, and transverse mass m^, where 
mj_ = \Jp\+ iti 2 for a particle of mass m. Defining ra- 
pidities and transverse momenta relative to the z axis, a 
4- momentum p^ = (E,p) can be written in two equiva- 
lent forms as 



(mj_ cosh y,p±, m± sinh y) 



m2 + P± cosh rj,pj_,p± sinh 77 



(4) 



In terms of the polar angle 6 from the z axis, ?y = 
— lntan(0/2). The standard velocity of a relativistic par- 
ticle is v = \p \/E. In our analysis, an important role will 
be played by the "transverse velocity" r defined by 



P± 



V- 



(5) 



In general, the pseudo-rapidity and velocity of a particle 
can be expressed in terms of r and y: 



v = v( r , y) = in 



y/ r 2 + sinh 2 y + sinh y 



(6) 



v(r,y) = 



\J r 2 + sinh 2 y 
cosh y 



For massless particles, r = v = 1 and r\ = y. 

Our focus will be on dijet event shapes in e + e~ colli- 
sions, which have the property that e — > implies back- 
to-back pencil-like jets (and e = for the lowest order 
partonic configuration e + e~ —> qq). To maintain simplic- 
ity, we will not include recoil sensitive observables such 
as broadening (67j in our analysis^ The event shapes 
we consider are also bounded as < e < e max , where 
the maximum value e max depends on the observable in 
question but is typically 0{\). For the e -C e max limit, 
one can sometimes simplify the expression defining e by 
neglecting corrections of 0(e 2 ). We will distinguish ex- 
pressions for event shapes that are valid under this ap- 
proximation by adding a bar, e. 

Various examples of dijet event shapes are shown in 
Tab. HI including their original definitions e and expres- 
sions e valid when e <C e max . For angularities and 2- 
jettiness, there are no simplifications in the dijet limit, 
so e = e without higher order terms. The notation for 
various items in the table require explanation. The nor- 
malization factors are defined by 



Qp = \& 



Q = J2 E > 



(7) 



2 The original thrust variable is T = 1 — 

3 The shape parameter H2 introduced in [6 
the C-parameter with the substitution Q 



661 is equivalent to 



4 For a recoil sensitive observable, the axis used to compute e 
(typically the thrust axis t) differs from the axis which minimizes 
e by an amount that can have an 0(1) effect on the value of e. 
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/ e (r, y) 


T 




T (a) 


c 


P± 


Original 


re -M 


e -\y\ 


L e -\v\(l-a) 
V 


3r 
cosh. Tj 


0(±y)e™ 


P-scheme 


re -M 


re -M 


re -M(l-a) 


3r 
cosh 77 


rO(±ri) e Tv 


E-scheme 


L e -M 

V 


L e -M 

V 


L e -i»jia-o) 

V 


r 3 
v cosh 77 


-6»(±r ? )e T " 

V 


R-scheme 


re -\y\ 


re -\y\ 


re -|»IU-<0 


3r 
cosh y 


r6{±y)e* y 


J-scheme 


e -\y\ 


e -\y\ 


e -|wl(l-o) 


3 

cosh y 


8(±y)e*y 



TABLE II. The functions f e (r,y) of "transverse velocity" r and rapidity y for various dijet event shapes using several different 
schemes for treating hadron mass effects. Here r\ = rj(r,y) and v = v(r,y) are given in Eq. ((6J. 



The unit vector t obtained in the minimization defining 
thrust r is referred to as the thrust axis, and we have 
taken it to be aligned with z to define rapidities. The 
angle between particles i and j is denoted by 9ij , and the 
angle between particle i and t is denoted by Qi. Finally, 
the ± labels refer to each of the two hemispheres defined 
by the plane normal to the thrust axis. For later conve- 
nience, we have written all e formulae in terms of p± , m± , 
y, and/or rj. Note that we can always replace Q p — > Q 
in the overall normalization for e, since the correction in 
doing so is beyond the order we are working. 

For a dijet event shape e, the largest part of the cross 
section comes from e -C e max , as do the most important 
hadronization corrections which are the focus of this pa- 
per. The simplicity of the dijet limit makes it possible 
to derive factorization theorems^ for these cross sections 
which facilitate calculations of higher order perturbativc 
corrections oc a^(ln fc e)/e, as well as defining nonpertur- 
bativc corrections in terms of field-theoretic matrix ele- 
ments. For our purposes, the relevant point is that at 
leading order in the nonperturbative corrections, we can 
split e into pcrturbative (e p ) and non-perturbative (ca) 
contribution^ 



e A 



(8) 



The e p term is generated by particles with momenta 
p Aqcd and here we can neglect corrections from 
hadron masses m# up to second order in the mn/Q 
expansion. The e\ term involves corrections from soft 
particles with momenta ~ Aqcd, and from Tab. HI we 
see that e\ ~ A/Q where A ~ ran ~ Aqcd- Thus simple 



power counting dictates that to determine e\, we cannot 
neglect hadron mass effects in the definition of e. 

For our analysis, we will find it convenient to charac- 
terize each event shape by a function f e (r, y) of the trans- 
verse velocity r and rapidity y, which we define from the 
e <C e max limit via 



e = Q^2 m i fe(n,yi) 



(9) 



The various examples in Tab.|T]have the following f e (r, y) 
functions: 



fr(r,y) = \Jr 2 +sinh 2 (y) -sinh|y| 
U 2 {r,y) = e~\y\, 

U w (r, y) = r a fvWshm 2 / _ 



(10) 



fc(r,y) 



dy 
3r 2 



l-a 



\/ r 1 + sinh 2 y 



(11) 



f p± {r,y) = 6{±y)e^ , 
where 

E coshy drj(r,y) 1 

\p\ yjr 2 + sinh 2 y dy v ' 

The translation from the notation of Ref. [l[ (de- 
noted with superscripts SW) to our notation is 
f! w (y,my P l) = f e (r,y)/r. 

B. Hadron Mass Schemes 



5 The simplest examples of dijet factorization rely on being able to 
write the dijet event shape as a sum of contributions from ener- 
getic collinear particles in the ± hemispheres (n,n), soft pcrtur- 
bative particles (s), and soft nonperturbative particles (A), via 
e = e n + e„ + e s + e a • See for example Ref. 16311 . 

6 The heavy jet mass event shape pjj = niax{p+,p_} does not 
admit such a decomposition, and correspondingly its factoriza- 
tion formula is a bit more complicated. The light jet mass 
PL = min{p+,p_} is not a true dijet event shape since pz, — > 
does not imply a dijet configuration. 



In experimental analyses, different "schemes" are often 
adopted for the treatment of hadron masses depending 
on the available information. In the context of event 
shapes, a detailed discussion of these schemes is given in 
Ref. [lj. These schemes correspond to different choices 
for f e (r,y) that yield the same value for the perturba- 
tive contribution e p but potentially different values for 
the nonperturbative contribution e^- The schemes con- 
sidered in this paper are summarized in Tab. [TXJ 
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In the "P-scheme" , one makes measurements with only 
3-momentum information. One performs the substitu- 
tion Ei — > \pi\ in the formula for event shapes, and cor- 
respondingly m± — > r mi in Eq. ^ . To satisfy infrared 
safety, the original f e {r, y) must tend to a constant value 
(possibly zero) in the r — > limit, and this implies that 
the P-scheme event shapes will always vanish linearly 
with r. The P-scheme replacement affects the jet masses, 
angularities, and 2-jettiness which become 



W P i 

<) = 7T E 1*1 ( sin ^) a (l - I cos^l) 1 -" , 



(12) 



P± = Ql 



E \a 



2 w i] 



(i,i)e± 



Another scheme discussed in Ref. [l| is the "E-scheme" 
(see Refs. [H |H [H for earlier discussions), where 
only measurements of energies and angles are used to 
construct observables. Compared to the P-scheme, one 
makes the substitution 



- ^ Ei _ 

Pi -> T-7T Pi 
\Pi\ 



(13) 



in the formula for event shapes. This modifies all exam- 
ples in Tab. U except for angularities: 



T = T 9 



I ^£,(1-1 COS 0*1): 



(14) 



EiEj sh 



2 "ij 



i,3 



2 ' 



sin 



Note that thrust and 2-jettiness arc identical in the 
E-scheme, = t e . In the P-scheme, the event shapes 
are all linear in r, and hence the corresponding E-scheme 
results for / e (r, y) are simply obtained by multiplying by 
E/\p \ = 1/v. Note that the E-scheme and P-scheme are 
defined in terms of pseudo-rapidity r\ (equivalently, the 
polar angle 9). 

In order to consider a wider range of observables, we 
will introduce two new schemes. In the "R-schcmc" 
(rapidity scheme), we take event shapes defined in the 
P-scheme and make the replacement r] — > y. To de- 
fine R-scheme event shapes where e ^ e, we carry out 
this replacement for e, and then define e R = e R . The 
"J-scheme" is the closest to the jet mass observables, 
and is defined by taking the R-scheme result and setting 
r = l, f e -'(r,y) = f e n(l,y). 

We emphasize that the naming of schemes discussed 
here is set simply by convention. For understanding 
power corrections one only needs to know the functional 
form of f e (r, y). 



C. Effect of Power Corrections on Cross Sections 

For recoil- less dijet event shapes that satisfy Eq. ((5J, 
the perturbativc/nonperturbative factorization of the dif- 
ferential distribution in the e — > limit implies 

£ = /<( e -|)^+<^ ,15) 

Here da/de is the most singular perturbative cross sec- 
tion to all orders in a s , and contains the full leading 
power perturbative soft function. F e is the shape func- 
tion that depends on the specific event shape one is in- 
terested in. It contains nonperturbative power correc- 
tions (and, as we will see, perturbative corrections). If 
Qe ~ Aqcdj then the entire function F e {l) has an im- 
portant impact on the cross section and in practice one 
models it with a few coefficients which can be fit to 
data [13, IHI , or uses a complete basis which can be sys- 
tematically improved (58j . 

For Qe ^> Aqcd, the function F e (£) can be expanded 
for t 3> Aqcd in terms of nonperturbative matrix ele- 
ments of operators. The first terms are 



F e (£) = 5(£)-8'{l) 9.\ + 



a s A 



QCD\ 

e ) 



(16) 



where fif (//) is a dimension- 1 nonperturbative matrix el- 
ement (defined here in the MS scheme) that encodes the 
power corrections we wish to study. It is defined by 



n; = (o | y^(Qg)y n y fl | o; 



(17) 



where (Qe) is a Q-indcpendent field-theoretic operator 
that measures the combination Qe, and Y (Y) are Wil- 
son lines with gluon fields in the fundamental (anti- 
fundamental) color representation along the directions 
specified by n = (1, t) and n = (1, —t). For example, 



Y n = P exp 



ig j dsn ■ A(ns) 
o 



(18) 



where P stands for path-ordering and is the gluon 
field. The fact that the dimension-1 matrix element f2f 
is nonperturbative is easy to understand on dimensional 
grounds since the only scale for this QCD vacuum matrix 
element is Aqcd- In Eq. (fT6|) the C(q; s Aqcd/^) term in- 
volves perturbative corrections to the leading power cor- 
rection which will be discussed in Sec. IV Bl 

Plugging Eq. (|TB)) into Eq. (fTS)) for the event shape 
distribution one finds 



da 
de~ 



da 
d7 



Of d da 
Q de de 



(19) 



where the ellipses denotes higher order terms in a s and 
Aqcd/^- Eq- (H3) corresponds to a shift e — > e — Of /Q to 
first order in l/Q, and reproduces the known shift found 
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in the dispersive approach [31H33I ]. Thus, the dominant 
effect of power corrections (hadronization) on dijet event 
shapes for Qe Aqcd is simply a shift in the distribu- 
tion. 

Following Ref. [3(| we note that as long as d<r/de tends 
to zero in the far tail region, one can derive an operator 
product expansion for the first moment of the distribu- 
tion as well. Defining the full and perturbative moments 



energy-momentum tensor T^ u (t,x) by [61M631 16 



Ida 



/ pert 



/ de e — — 

a de 



(20) 



one can use the factorization in Eq. (fT5|) and the expan- 
sion of the shape function in Eq. (fT()]) to show 



/pert 



~Q 



(21) 



where again the ellipses denote higher order terms. For 
the event shapes in Tab. U all the event shapes except 
for C tend to zero in the far tail region, so the leading 
power correction generates a Q-dcpendcnt shift of their 
first moment. We will use this feature to extract flf 
from the Pythia 8 and HerwigH — h hadronization models 
in Sec. IV D[ See Ref. 47J for a discussion of the modifi- 
cation necessary for the C-parameter moment. 



III. THE TRANSVERSE VELOCITY 
OPERATOR 

The goal of this paper is to study the effect of hadron 
masses on the power correction in Eq. (|17p . In order to 
formulate the operator (Qe) appearing in the definition of 
fif , we will follow the energy-momentum tensor approach 
of Rcfs. 0, H3|, and generalize it so that we can treat the 
dependence on the transverse velocity r. In particular, 
the event shape e will be written as the eigenvalue of an 
operator acting on the final state which includes hadron 
mass effects. In Sec. IIII B[ we consider the crucial role of 
boost invariance for identifying universality classes. 



A. Comparison to Transverse Energy Flow 

To set up our analysis, it will be convenient to first re- 
view definitions from the literature that do not account 
for hadron mass effects in the event shapes. These cor- 
respond to setting r = 1 for the event shapes defined by 
Eq. ([5]). The transverse energy flow operator £t{v) is 
defined by 



S T ( V )\X) = ^2pt5( V -r )i )\X) 



(22) 



i£X 



(Here and below we suppress the dependence on the 
thrust axis t.) The operator £t(v) is related to the 



1 



cosh r/ 



lim R 2 dth t T oi (t,Rh), (23) 



where n is a unit vector pointing in the (8, (f) direction 
for the 8 corresponding to rj^\ Using £t(?7), we can define 
an operator for the event shape e by Q 



e Q \X) 



1 

Q 



d v f e ( v )£r(v)\X) 



(24) 



where in the notation defined in Eq. @, this 
fe(v) = fe( r = U V — v)- F° r massless particles with 
r = 1, this operator satisfies io\X) = e(X)\X), so the 
measurement operator in Eq. (|17|) is given by Ai e = Q&o- 
To generalize this operator formalism to include 
hadron masses, we require a transverse momentum flow 
operator that is more differential, namely a y) 
which can pick out states that have a particular trans- 
verse velocity r. We will refer to this simply as the "trans- 
verse velocity operator" , and define it by its action on a 
state 

£ T (r,y)\X) = Y,mi5(r-r i )5{y-y i )\X) . (25) 

iex 

Note that here we use rapidity y rather than pseudo- 
rapidity rj. In App. |X]we show that £r(r,y) can be de- 
fined in terms of the energy-momentum tensor as 



£r(r,y) = 



sech y 



yj r 2 + sinh 2 ?/ 



lim R 

R— > 00 



2tt 



hi T (R,Rvh) , 
(26) 



where v = v{r,y) and r\ — T}(r,y) are given in Eq. ([5]). 
The unit vector h again points in the (8, 4>) direction and 
hence depends on y and r through its dependence on 
7] = -lntan(0/2). 

The physical picture for the distinction between £t(v) 
and £t{t, y) is shown in Fig. [T] The energy flow opera- 
tor £t(v) involves an expanding sphere of radius R inte- 
grated over all time, and measures the total transverse 
momentum for rapidities in an infinitesimal interval Srj 
about rj. The transverse velocity operator £t{t, y) in- 
volves a spheroid that expands in both space and time 
with a finite velocity v, and it measures the total trans- 
verse mass for particles in an infinitesimal interval in both 
77 and the velocity v (or equivalently an infinitesimal in- 
terval in y and r). Using £t{t, y), the value of an event 



In the proof of Ref. [13 that Eq. (J23J reproduces Eq. {22j for 
scalars and fermions, one assumes that all particles are massless. 
If the energy momentum tensor is considered for massive fields, 
then Eq. II23H yields Eq. 1221 1 . but rj must be identified with 
pseudo-rapidity and the factor pr- — > mj". 
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y) 



FIG. 1. The energy flow operator (top) compared to the trans- 
verse velocity operator (bottom). Measurements are made 
with respect to the thrust axis t. The dashed arrows corre- 
spond to particles with lengths given by the particle velocities. 
Shading indicates which particles are measured by the opera- 
tor. Note that the velocity v(r, y) and pseudo-rapidity r)(r, y) 
are functions of the transverse velocity r and rapidity y. 



shape e for a state \X) with massive or massless particles 
is given by e\X) = e(X)\X) where the operator 



1 



+ 00 



dy / dr f e (r,y)£ T {r,y) 



(27) 



involves f e (r, y) defined in Eq. (|9]) . This is the desired 
generalization of Eq. (|M|) that will allow us to treat the 
effect of hadron masses on event shape power corrections. 
The result in Eq. (|2"T1) completes the matrix element def- 
inition of Of given in Eq. (fT7|) . 



Due to the invariance of the vacuum |0) and Wilson lines 
Y n and Y ri under this boost, we can choose y' = —y so 



(0 | Y{Y^£ T {y)Y n Y n | 0) = (0 | YlY^£ T (0)Y n Y n | 0). 

(29) 

Using Eq. (|24[) . again with r] = y, the power correction 
in Eq. (|17|) simplifies to 



+ 00 



dvfeiv) 



Q\Y{Y^£ T (0)Y n Y n \0) 



(30) 



Thus as argued in Ref. Q , the power correction is univer- 
sal if one neglects hadron mass effects, since it depends on 
a common nonpcrturbative matrix element times a cal- 
culable r\ integral specific to an event shape. This result 
agrees with the the dispersive approach [3ll - l33| , and also 
explains why it only applies for the first power correction 
in the expansion of Eq. (jT6j) . 

We can apply the same logic to our transverse velocity 
operator in Eq. (f25j) . now accounting for the effect of 
hadron masses. Under a boost of rapidity y' , 



U(y')£ T (r,y)U(y'^ = £ T (r,y + y') 



(31) 



Choosing y' = —y, we find that the leading power cor- 
rections to dijet event shapes are all described by the 
nonpcrturbative matrix element 



fli(r) = {0\Ylrtir(r,0)Y n Yn\0) 



(32) 



which depends only on r and is independent of the event 
shape e. Using Eq. (f2"T)) . the power correction Q,\ simpli- 
fies to 



dr 



dyf e (r, y) 



fii(r). 



(33) 



We will consider the implications of Eqs. (|32p and (|3"3")l 
for universality in the next section. Note that there is no 
limit where Eq. (|33|) approaches the massless approxima- 
tion in Eq. (|3"DJ). 



IV. UNIVERSALITY FOR EVENT SHAPES 



B. Boost Invariance 

Both £t(v) an d ^T( r -,y) have nice transformation 
properties under longitudinal boosts. These arguments 
were first given in Ref. in the context of £t(ij) to prove 
universality of power corrections for massless particles, 
and we will extend the logic for £x(r,y) to develop the 
notion of universality classes which account for hadron 
masses in the next section. 

For the case of massless particles in Ref. , f] — y and 
it was shown that under a boost of rapidity y' along the 
thrust-axis: 



U{y')£ T (y)U{y')^£ T (y + y'). 



(28) 



A. Universality Classes Denned by g(r) 

The boost invariance logic of the previous section 
shows that for the power correction VL\ , we can factor out 
the rapidity dependence from the nonpcrturbative matrix 
element. We define the integral appearing in Eq. (fJ3"|) as 



(34) 



dy fe(r,y) = c e g e (r) , 
where ff e (l) = 1 and 

/+oo 2 r+°° 

dyf e (l,y), 9e(r) = — dyf e (r,y). (35) 
-oo J — oo 
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T 


T-> 


T(o) 


c 




Common 


2 


2 


2 

1-a 


3tt 


1 



TABLE III. Expression for the c e coefficients for various dijet event shapes from Table [TT1 Since c e are defined using fe(l,y), 
they have the same value in each class. 



9e(r) 


T 






C 


P± 


Original 


9r{r) 


1 


C 


^X(l-r 2 ) 


1 


P-scheme 


9r{r) 


9r(r) 


9V (r) 


^K(l r 2 ) 


Sr(r) 


E-scheme 


r 


r 


r 


r 


r 


R-scheme 


r 


r 


r 


r 


r 


J-scheme 


1 


1 


1 


1 


1 



TABLE IV. The functions g e (r) of transverse velocity r for various dijet event shapes from Table HIl Event shapes with the 
same g(r) belong to the same universality class. 



(For the special case of event shapes that vanish for mass- 
less partons, J dyf e (l,y) = 0, we define c e = 1 and let 
<7e(l) = 0.) Thus, the leading power correction for an 
event shape e can be written as 

Q? = c e Of, Of = [ drg^n^r), (36) 
Jo 

where Oi(r) is given in Eq. (|32l) . 

Using the fact that for massless hadrons f e (v) = fe( r = 
l,y = j]) in Eq. (|3T)1) . we see that the coefficients c e are 
precisely the classic universality prcfactors that one de- 
rives neglecting the hadron mass dependence of event 
shapes 0, H HI HI HI El H-fzl- The function g e (r) 
then encodes the effect of hadron masses through the 
nonperturbativc parameter Of . 

The key result from Eq. (|36p is that each unique func- 
tion g e (r) defines a universality class for dijet events 
shapes. In particular, for two different event shape vari- 
ables a and b, if g a {f) = .9b(?') then Of = Of (equiva- 
lently f2°/c a = Cl\/cb), so their power corrections agree 
up to the calculable constants c a and Cb- We say that two 
such event shapes belong to the same universality class, 
and we will often refer to Of as the universal power cor- 
rection defined by g{r). 

Recall from Tab. Q] that the f e ( r ,y) functions in gen- 
eral depend on the measurement scheme used for treat- 
ing hadron mass effects (E-scheme, P-scheme, etc). The 
functions g e (f) will also in general vary with the mea- 
surement scheme. However, the c e coefficients arc inde- 
pendent of the scheme for treating hadron masses since 
they are defined with r = 1 corresponding to the mass- 
less limit, and f e (r = 1, y) is the same in all schemes. Our 
classification of universality classes with common g(r)'s 
is the same as the identification of event shapes that have 
the same hadron mass effects made in Ref. and we 



will elaborate on the precise notational relationship in 
Sec. HVC] below. 

A summary of coefficients c e , functions g e {r), and uni- 
versality classes are collected in Tables IIII1 IIV1 and |V] 
and will be discussed in detail in Sees. IIVBI and IIVCI 
below. Since many of the standard event shapes have 
different g(r) functions, their power corrections are not 
related by universality. We will take up the question 
of numerically approximate relations between power cor- 
rections in different universality classes in Sec. IIVDI Fi- 
nally in Sec. IIV El we consider the impact of changing 
the renormalization scheme defining Oi(r) from MS to a 
renormalon-frce scheme. 

B. Generalized Angularities 

In order to see how universality works in practice, it is 
instructive to consider a family of event shapes which are 
a simple generalization of angularities, and are labeled by 
two numbers n > and a < 1: 

i 

This corresponds to extending the R-scheme definition 
of angularities by incorporating a positive power n of 
r = pj_/m±, and this n dependence allows the event 
shape to directly probe hadron mass effects. In terms of 
Eq. ([9]) one trivially finds 

f n ,a(r,y) = r n e-^ 1 - a K (38) 

These generalized angularities all have the same value of 

/ + OO Q 
dye -\v\(i-a) = -— = Ca , (39) 
-oo -L a 
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which from Tabic Mil is also the same as for classic angu- 
larities T( a ) . Computing the function encoding the mass 
dependence we have 

(l-o) 



9n,a{r) = 



dy r n e 



(40) 



Thus the Ti n ,a) event shapes belong to universality classes 
labeled by n, and represented by the functions 



9n(r) 



(41) 



Each value of n defines a different universality class, so 
in general there are infinitely many different event shape 
universality classes. For the r™ class, the universal power 
correction from Eq. (|36l) is 



drr n ili(r). 



(42) 



This procedure, of generalizing an event shape to ob- 
tain different sensitivity to hadron mass effects by multi- 
plying by r n , can just as easily be applied to event shapes 
other than angularities. By choosing large values of n, we 
preferentially select out particles whose transverse mo- 
mentum is large compared to its mass. This might be 
useful in an experimental context to deweight the contri- 
bution of soft particles to event shapes, beyond the linear 
suppression in m± necessary for infrared safety. 

We have included a summary of some common uni- 
versality classes in Table fVl For n = and n = 1, the 
universality classes correspond with two that appear for 
classic event shapes, so we will also refer to g{r) = 1 as 
the Jet Mass class and g(r) = r as the E-scheme class. 

In Sec. lVDl we will show that Pythia 8 and Herwig++ 
do indeed exhibit universality when holding n fixed but 
varying a, but as expected yield different power correc- 
tion parameters O™ when n is varied. 



C. Classic Event Shapes and Mass Schemes 

In this subsection, we discuss c e and g e (f) for the 
more traditional event shapes enumerated in Tab. HI and 
show how g e (r) changes when using various measurement 
schemes for hadron masses. Results for the correspond- 
ing f e (r,y) were summarized above in Table [TT1 As al- 
ready mentioned, the results for c e are independent of 
the measurement scheme for hadron masses, and can be 
computed directly with Eq. ([55"]). For the various classic 
event shapes they are summarized in Table Hill 

For the event shapes in Tab. Q] with their original defi- 
nitions, integrating their f e (r,y) over y we find 



9r(r) = 


1 - E(l 


: 9p± 0) = 


1, 


9r ia) (r) = 


r, 


9c(r) = 


2r 2 


—K(l 

7T 



where E and K are the complete elliptic integrals 

fir/2 



pit/* 

E(x) = / d0(l-a;siir 
Jo 



2^1/2 



(44) 



it/2 



K(x) = / d0(l-xsiiT 0) 
Jo 



2/U-1/2 



A plot of the g{r)'s in Eq. 03]) is displayed in Fig. [2] 

Direct analogs of _g T (r), gc{ r ), 9p±{ r ) were computed 
in Ref. [l[. The translation from their notation (super- 
scripts SW) to ours is cf w = c e and 5cf w (m 2 /p 2 ± ) = 
c e [g e (r)/r - 1]. Fig. 1 of Ref. plots (1 + 5cf w /cf w ) 
versus pj_/m, which is the direct analog of our Fig. [2] 
In our figure, the function is bounded because of our use 
of the r = p±/m_i_ variable rather than pj_/m. These 
bounded g(r) functions are more convenient for our ba- 
sis discussion in Sec. IIVD] below. 

By looking for event shapes with common g(r)'s in 
Eq. (|43|) or Fig. [2] we see that p± and t 2 are in the "Jet 
Mass" universality class, and that the angularities tv ) for 
any a belong to a "E-scheme" universality class. Defining 
ill an d H\ as the universal power corrections for the Jet 
Mass and E-scheme classes, respectively, and accounting 
for the c P factors, we have 



ill 2 =20?, 



ny = 



1 -a 



Of 



(45) 



Thrust and the C-parameter have their own g(r) func- 
tions and hence among these event shapes are alone in 
their universality classes, with power corrections il\ and 

We now consider the various hadron mass schemes. 
For all of the event shapes in Tab. U in the E-scheme, 
one has a common g e {f) = r, so they all belong to the 
"E-scheme" universality class. It was for this reason that 
Ref. [l[ considered the E-scheme to be privileged hadron 
mass scheme. In this sense, universality in the E-scheme 
is closest to the universality for massless particles, since 
any event shape that can written as a function of four- 
vectors has an E-scheme definition and they all have the 
same universal power correction. For event shapes de- 
fined directly in terms of m±, r, and y (e.g. the gen- 
eralized angularities in Eq. ([37])), the E-scheme can be 
defined by first expressing the event shape in terms of 
four- vectors (with no explicit m dependence), and then 
applying the E-scheme replacement in Eq. (|13|) . Then 
any E-scheme observable has g(r) = r and hence is in 
the E-scheme class. An example of an exception are event 
shapes that vanish in the massless limit, g e (r = 1) = 0, 
which do not have a meaningful E-scheme definition. 

From Table HVl we see that the event shapes mea- 
sured in the R-scheme also fall into the E-schcmc class 



8 Accounting for the c e dependence, the universal power correc- 
tions for the thrust and C-paramctcr classes arc Q® T = f2J/2 
and nf c = nf/3n\ 
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Class 




Event shape 


Jet Mass class (fi? or Q J) 
E-scheme class or fif ) 
r" class (fi?) 


1 

r 

r° 


P±, t"2, r , r (a) , C J 

T(a), T —T 2 ,C , p±, T , T (a) , O , p ± 

generalized angularities r(„ ja ) in Eq. (|37p 


Thrust class (fi? T ) 
C-parameter class (fif c ) 
r 2 class (OJ) 


9oW 
r 2 


p p 

c 

p 

T (2,a), T( a _ ) ._ 00 ) 



TABLE V. Event shape classes with a universal first power correction parameter ilf" . For a given event shape, the full power 
correction is f2f = c e ilf° . 



(g(r) = r), while event shapes measured in the J-scheme 
fall into the Jet Mass class (g(r) = 1). The fact that 
E-scheme and R-scheme jet shapes are in a common uni- 
versality class is a non-trivial consequence of the relation 
dij/dy = l/v from Eq. (p~T|) . 

For the P-scheme, universality classes are more com- 
plicated. Carrying out the integral in Eq. ([55)1 to find 
the g(r)'s one finds that thrust and the C-parameter 
are unchanged. For the jet masses and 2-jettiness 
g p p (r) = g T p (r) = g T (r) , so both p± and belong to 
the thrust universality class: 



1 



1 



1 o 1 r) 1 I 



(46) 



For P-schcme angularities, g T p(r) docs not appear to 

have a simple analytic formula for arbitrary a (although 
it is easy to compute numerically). By making a change 
of variables y = — lntan(x/2), a convenient way of writ- 
ing it is 



( r ) = (!-«) / d X 
{a > Jo 



1 — (1 — r 2 ) sin x + cosx 



r a 2 sin a x 



a-l 



(47) 



For integer values, it is simple to find an analytic form, 
and for a = — 1, —2 we finq^l 

g T P_ i} (r) = 2 - r - 2 Q- - r) ln(l + r) , (48) 



9r r _ 2) (r)=9-^-{7-^)E(l-r 2 ) 

- (4-3r 2 )A"(l-r 2 ). 

For any a, g T p (0) = and g T p^(l) = 1. For large nega- 
tive values of a one finds 



(49) 



9 In general for odd a, g(r) involves only ln(l + r) and for even a 
only the elliptic functions K(l — r 2 ) and E(l — r 2 ) appear. 




FIG. 2. From top to bottom the g e (r) functions for Jet Masses 
and 2-jettiness, Angularities, C-parameter, Thrust, and the 
P-scheme angularity with a = — 1. Lines correspond to exact 
numerical values, open circles to two terms from the basis of 
Sec. IIVDI and filled circles to three terms from the basis. 



Hence angularities in the P-scheme for large negative a 
belong to the same class as the generalized angularities 



in Eq. fl37|) for 



2, with fi/' 



T (2,a) 



In practice, g T p^(r) for arbitrary a quickly converges 
towards g T p (r), and hence there is a quasi-universality 

(— oo) 

for angularities in the P-scheme. Also gc (r) and g T (r) are 
not so different from g T p ^ (r), implying an approximate 

universality between event shapes in different classes for 
all event shapes in the P-scheme. This was already noted 
for thrust and the C-parameter in Ref. [l|. We will next 
develop a complete basis for describing g e (r) functions 
that will allow us to make this observation more quanti- 
tative. 
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L 

00 


01 


02 


03 


Jet masses, T2 


1 











Angularities 


1 

2 


1 

2~7I 








Thrust 


0.383 


0.299 


0.050 


-0.006 


p 
T (-D 


0.355 


0.295 


0.064 


—0.004 


C-parameter 


0.393 


0.300 


0.046 


-0.007 


~P 


1 


1 


1 





' (a— > — oo) 


3 


2^3 


6\/5 



TABLE VI. Numerical value of the coefficients of the complete 
basis for the various g e (r) functions. 

D. Orthogonal Basis for fii(r) 

Many g e (r) curves are still parametrically close, even if 
the corresponding event shapes are formally in different 
universality classes. Examples are thrust, C-parameter, 
and f^x) angularity shown in Fig. [2] To get a quanti- 
tative handle on this observation, we can use a complete 
set of orthonormal functions for r G [0, 1]: 

h n (r)=V2^+TP„(2r-l), (50) 

/ drh n (r) h m (r) = 6 nm , 
Jo 

where P n are the Legendre polynomials. Now we can 
decompose any of the g e (r) functions in this basis in the 
usual way: 

oo 

9e(r) = jT l b°h n (r), (51) 

n=0 

K = / drg e (r)h n (r) . 
Jo 

Since the g(r) functions for classic event shapes are fairly 
close to low-order polynomials, the first few terms in the 
basis will provide an accurate approximation. 

The values of the b e n coefficients are shown in Table IVII 
for the classic event shapes. The approximation for the 
g e (r) functions are plotted in Fig. [5J where the exact 
values are shown by lines and the approximation with 
two terms from the basis (bo, b±) are shown by hollow 
circles, and with three terms (bo, &i, b 2 ) by filled cir- 
cles. For the jet masses, 2-jettiness, and angularities, the 
approximate result is exact with two terms in the basis. 
For the remaining events shapes, the approximation with 
two terms is likely sufficiently accurate at the level one 
expects of current experimental and perturbative preci- 
sion. With three terms, the approximation is excellent 
in all cases, so the third term can be regarded as a high- 
precision correction. This is also apparent from the val- 
ues in Table IVII where bo and b\ are much larger than bi 
(and computing b n> 2 one finds they are negligible). 



Using Eq. ([50]) one can write any Of in terms of a 
denumerable set of power correction parameters: 

oo -i 

Slt = J2Kn[ n) , n[ n) = drMr)fii(r). (52) 

Only the first few terms will be numerically relevant for 
most event shape observables. Using Eqs. (f5"Tj) and (f52|) 
and the results in Table IVII we can make the following 
exact identifications 

= of , a* 1 * = 2V3 of - V5 of . (53) 

Thus the power correction parameters for the Jet Mass 
class and E-Scheme class already give an excellent ap- 
proximation for the various classic event shapes. To re- 
fine the prediction even further one can use the next term 
in the basis, 0^0 Writing the leading power correction 

for the other event shapes in terms of Of, Of, and O^ 2 ' 
we have: 

Cll = 1.034 Of - 0.135 Of + 0.050 f4 2) , (54) 
Of = 1.039 Of - 0.127 Of + 0.046 0^, 

p 

Oi C_1) = 1.022 Of -0.156 Of + 0.064 O^ 2 } . 

In general both the Of and Of terms are numerically 
important since experiment favors values Of ~ Of/2, 
(2) 

and O^ ; can be typically neglected. However, since the 
numerical coefficients in Eq. (fM)) are so close, one is jus- 

tified in using the approximation 0[ ~ Of ~ fl 1 1 up to 
corrections of ~ 15%. 

We note that the above basis analysis is valid only 
at a fixed value of Q since Oi(r, fi) has an anomalous 
dimension which we will compute in Sec. [Vj and the ap- 
propriate n scales with Q. As we will see in Sec. IV C| we 
will need to refine the above analysis in order to relate 
power corrections at different values of Q. 

E. Renormalon-free Definition of fl±(r) 

It is well known that the first moment of a perturbative 
event shape distribution has a Aqcd renormalon ambigu- 
ity in the MS scheme [72| , which corresponds to a renor- 
malon in the MS power correction parameter Of that we 
have been considering so far. To our knowledge, the only 
renormalon-based analysis of event shapes that observes 
sensitivity to hadron masses is Rcf. [34[ (which in turn, 
called into question the massless universality results). In 
this section we argue that the universality relations given 



In principle, we could extract exactly from the r ^_ > _ 00 ) 

power correction, since from Table IVTl we see that g p (r) 
is saturated by the first three terms in the Legendre expansion. 
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in Sec. lIV Al will remain unchanged as long as one defines 
appropriate renormalon-free schemes for the 1/Q power 
corrections. 

In general the Aqcd renormalon is removed by con- 
verting the power correction parameter to a new scheme 



(55) 



where S e (R,(j>) is a series in a s ([i). There is a correspond- 
ing change to the perturbative part of the cross section 
that depends on this same series. Writing a e (x) for the 
Fourier transform of da/de this change is 



(56) 



Recall that universality classes relate Of (fi) for different 
event shapes e, and that the relations are nonperturba- 
tive. Hence it is clear that one has the same renormalon 
ambiguity for event shapes that are members of the same 
class (differing only by the proportionality constants c e ). 
Thus it is sufficient to adopt a common scheme change 
for members of the same class via any scheme satisfying 



6 e (R,n)=c e Re' 1E 



E 

n=l 



(^))>.) Wfl) . (57) 



Here the 5n°'\n/ R) coefficients involve factors of 
ln(/i/i?) and are chosen such that the new £1\{R, is 
free of the leading Aqcd renormalon ambiguity. The nu- 
merical values for 5n'' ^ will depend on the representative 
e' that we choose, but any representative has the same 
renormalon and hence is a valid subtraction series for all 
members of the g e class. 

One popular scheme for fixing the coefficients 

Sn £ '\n/R) is based on the dispersive approach, where 
R = hi [31M33I ]. In this framework to obtain an appro- 
priate subtraction at 0{a 2 s ) one must also include the 
Milan factor [3^, H3|. As as typical of scheme changes 
based on QCD perturbation theory, this scheme removes 
the renormalon for the massless limit r = 1. 

At higher orders in a s a more convenient scheme uses 
the leading power perturbative soft function itself [57} , 
S'P ort , which does not require additional computations 
beyond those needed for a resummed analysis of the event 
shape. Again this scheme change removes the renormalon 
for the massless limit r = 1 . For each g e universality class 
we pick a representative e' that is a member of the class. 
Then for all event shapes e £ g e we defined 



S e (R,fi) 



■Re lE 



d ]n(ix) 



\nS^ Tt (x,n) 



(iRe->E)- 



(58) 



11 To facilitate a multiplicative renormalization structure and to 
account for non-Abelian exponentiation, one uses the logarithmic 
derivative of the position space soft function S| lert (a;, fi) l73f . a 
definition adapted from the top jet-mass definition of Ref. 1741 . 



For e = t Eq. (|58| yields the subtractions for thrust de- 
fined in Ref. [73J. Given that renormalons probe the in- 
frared structure of amplitudes, one might naively expect 
that different subtractions would be necessary for the dif- 
ferent universality classes, since they treat hadron masses 
differently. However, the scheme change in Eq. (|58|) is 
based solely on the QCD perturbative calculation of S 1 ^,™' 
with r = 1, so only the c e coefficients will differ in sub- 
tractions for different event shapes. Thus in this setup, 
the same subtraction 5 e > (i?, fj>) can be chosen for all uni- 
versality classes. 

In App. [B] wc carry out a standard renormalon bub- 
ble sum calculation for an arbitrary dijet event shape e 
satisfying Eq. (j9]) and demonstrate that Eq. ([58]) yields 
a perturbative cross section a e (x) and power correction 
parameter Jli(i?, (J.) that are free from the Aqcd renor- 
malon probed by this method. It might be interesting to 
consider extensions of Eq. (f5"5|) that satisfy Eq. (|57p but 
have additional dependence on hadron mass effects. 



V. ANOMALOUS DIMENSION OF fti(r) 
A. Running at One-Loop 

In this section, we compute the one-loop anomalous 
dimension of the QCD matrix element f2i(r, fi) in the MS 
scheme, with details given in App.0 We regulate the UV 
with dimensional regularization (DR) using d = 4 — 2e. 

Since £li(r) has mass-dimension one and is propor- 
tional to the infrared (IR) scale Aqcd, one must be care- 
ful to establish IR regulators for the perturbative calcula- 
tion in such a way that there is nonzero overlap with the 
operator matrix element. This requires at least one di- 
mensionful IR regulator, as well as a mechanism to probe 
different values of r. With this, we can then compute the 
anomalous dimension just as we would for any external 
operator in QCD. The anomalous dimension will be in- 
dependent of the precise IR procedure used to identify 
the matrix element. 

A convenient choice for the IR regulator is obtained by 
coupling a massive adjoint background source Jt* A to the 
Wilson lines in f2i(r) by the replacement 



A" A (x) -> A^ A (x) + J^ A (x) 



(59) 



in Eq. (|T8|) and in the QCD Lagrangian. There is no 
Lagrangian mass term for the source J fJ " A , but for it to 
serve as an IR regulator, we will consider it to carry a 



massive particle momentum where q 



Recall 



from Eq. (f25|) that the Sr(r, 0) operator in £li(r) sums 
over contributions from individual particles in the final 
state. If £t{t, 0) acts on any particle other than J^ A 7 
then the corresponding phase space integral is scaleless 
and dimensionful, and hence vanishes in DR. Thus the 
first nonzero contribution occurs where £t{t, 0) acts on 
a J tiA (and we then set all other J^'s to zero). This 
IR regulator is convenient for bookkeeping and does not 
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FIG. 3. Tree level graphs for Qi(r,/i). 



overly complicate the evaluation of loop integrals. Ef- 
fectively it amounts to considering one of the final state 
gluons, namely the one acted upon by £t(t, 0), as having 
mass m and thus r ^ 1. For convenience when drawing 
Feynman diagrams, we use the same notation for mass- 
less gluons and the source, and simply note that we must 
sum over the cases where each final state gluon is the 
source. 

At tree level, as shown in Fig. [3l we have only the 
source line. This yields the nonzero matrix element 



M tree 



(r) 



2a s C F 



(1 -r 2 )! 



(60) 



Hence our setup provides nonzero overlap with the oper- 
ator in Qi(r). 

We now turn our attention to the 0(a s ) corrections. 
Here we must fix a gauge for the A^ A massless gluons. 
We have carried out all our calculations both with tradi- 
tional Feynman gauge, as well as with a background field 
Feynman gauge where the source J^ A takes the place 
of the external background field. Both yield the same 
results. In addition to J^ A , we must introduce extra 
IR regulators specific to individual diagrams. We find 
that shifting eikonal propagators involving the loop or 
phase-space momentum k^ by n ■ k — > n ■ k + A„ and 
n ■ k — >• n ■ k + A fi suffices to regulate other IR diver- 
gences. 

For the diagrams in Fig. U which involve two parti- 
cles in the final state through either a ghost bubble, or a 
gluon bubble, the sum of graphs does not involve a UV 
divergence (moreover these graphs arc IR finite). Hence 
they can be ignored for the purposes of calculating the 
anomalous dimension. The remaining diagrams are ei- 
ther Abelian with color factor Cf. or non-Abelian with 
color factor CfCa- Abelian contributions are shown in 
Figs. [5] and [51 and in the sum over Abelian diagrams, 
the real radiation and virtual contributions exactly can- 
cel (for all e). It is easy to prove that this cancellation 
happens to all orders in perturbation theory. This proof 
uses the fact that Y[A + J] = Y[A]Y[J] for an Abelian 
theory with no light quarks, and that yLAjy^LA] = 1. 

This leaves just the non-Abelian diagrams coming from 
Fig. |B1 and triple gluon vertex diagrams in Fig. [JJ There 
is also a contribution from gauge coupling renormaliza- 
tion which is not just canceled by the vacuum polariza- 
tion graphs in many gauges. For the 1/e poles, we can 
take A nt n — > in the sum of graphs in Fig. [5] and seper- 
ately in Fig. [TJ so the extra IR regulators cancel out of the 





FIG. 4. Diagrams involving gluon or ghost bubbles with two 
cut particles. The 4 additional diagrams obtained by a flip 
about the horizontal or vertical axis are not shown. Diagrams 
with one cut particle for wavefunction renormalization and 
coupling renormalization are also not displayed. 




FIG. 5. Purely Abelian 0(a^) diagrams. Either gluon line 
crossing the cut can be the source. The 4 additional diagrams 
obtained by a flip about the horizontal or vertical axis are not 
shown. 



UV terms as expected. The sum of diagrams in Fig. [5] 
and the sum in Fig. [7] each have 1/e 2 poles, but these 
cancel in the complete sum. This leaves only a nonzero 
1/e pole which will yield the anomalous dimension. 

The final result for the UV divergence in the matrix 
element of the bare fli(r) operator is 

M l-loo P(r) = ^ _ ^4 ln(1 _ ^ M trcc (r) _ (gl) 



We define the renormalized MS operator by 
VL\^{r) = Z{r,^e) Oi(r,/i), 



(62) 

so Eq. (|61"j) determines Z(r,fx,e) at 0(a s ). Using 
/ida s /d^i = — 2ea s + . . . , the one- loop anomalous di- 
mension of f2i(r) is 



/x-^-Oi(r,/i) 
d/j, 



ln(l 



fii(r,/i). (63) 



Note that this anomalous dimension is positive since 
ln(l - r 2 ) < 0. 

Intriguingly, the anomalous dimension is r-dependent, 
showing the important role of hadron masses. However, 
there is no mixing for operators at different values of r, 
so this RGE can be solved exactly to yield 



Qi(r,/j,) = Qi(r,/i ) 



a s {fJ,n) 



(64) 
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FIG. 6. Independent emission diagrams with Abelian and 
non-Abelian contributions. The 4 additional diagrams ob- 
tained by a horizontal flip or complex conjugation are not 
shown. 




FIG. 7. Triple gluon Y-diagrams for the t>(a 2 s ) correction to 
Qi(r). The 12 additional diagrams obtained by a horizontal 
flip or complex conjugation are not shown. Diagrams with all 
3 gluons coupled to Wilson lines of the same direction vanish. 



r> 1 (r, A/0 )[l-r 2 



2C_a 

3o 



s(^o) 



Here one can consider /io ~ 2 GeV as the low energy 
hadronic scale where we specify the nonperturbative ma- 
trix element, and fi as a high energy scale that is appro- 
priate for the observable being considered. In Scc. lVBl we 
will show that /i ~ Qe for the region of event shape distri- 
butions with e -C e max and Qe ^> Aqcd, and \i ~ Qe max 
for first moments of event shapes. 

In order to use the resummed expression for f2i(r, \x) 
to predict the evolution from Q| (/i ) to fif (/x), one would 
need to know the full r-dependence of f2i(r, /xo) to per- 
form the integral over g e (f)- We will see how to approx- 
imately circumvent this problem in Sec. IV CI 

One can also consider expanding Eq. (|64j) pcrturba- 
tively in a s (/xo) which yields 



fii(r,/x) = f2i(r,/x ) 



(65) 



In f JL) l n (i _ r 2 ) + 
V/i / 



If one truncates at 0(a s ), then one only needs two non- 
perturbative parameters defined at no to determine the 
the power corrections for an event shape at a higher scale 
xt: 



Qfifi) = SlXfto) + 



<^dl n fA) fi e,m (Mo) (66) 



Here fll(fio) is our standard power correction parameter 
at the scale /xo given by Eq. (|36p with Oi(r, /xo), and the 
slope parameter f2^' ln (/xo) is denned by 



^i' ln (Mo) 



-J dr ln(l 



, g e {r)Qi(r, /x ) . (67) 



The expanded form in Eq. (|66p is a reasonable approxi- 
mation if n is not too different from [1q. It works for a 
larger range than one would naively expect since there 
are numerical cancellations in the 0(a 2 ) term between 
ln 2 (l — r) and ln(l — r) contributions. 



B. The Wilson Coefficient of fii(r, /x) 

Having established that the nonperturbative matrix 
element f2i(r, fj,o) runs, we reconsider the operator ex- 
pansion of the shape function F e (£) in Eq. (|16p . now in- 
corporating a s corrections through a Wilson coefficient 
Cf(£, r, /j.) for f2i(r, n). The formula in Eq. (fTB]) becomes 



F e (£) = 6(£) + jdr Ct(e,r,n) c e5e (r) Ox^/x) 



(68) 



As usual the /x dependence of C\ (£, r, /i) cancels that of 
ili(r,fi). The dependence of Cf(£,r,fi) on £ and /x will 
determine the appropriate scale /x where there are no 
large logarithms in this Wilson coefficient. This in turn 
will determine the appropriate perturbative scale fi for 
the endpoint of the evolution derived in Eq. (|64[) . Since 
the £ dependence is treated differently by event shape 
distributions and by their first moments, a different scale 
/x will be found for these two observables. 

Taking Eq. (|?J3")) together with the cancellation of the 
/x dependence implies 



M A Cl e (€,r, M ) 



C A a s (iJ,) 



ln(l-r 2 )C^,r,/x). (69) 



At order a s using Eq. ([To]) this becomes 



^ T Ct{£,r,fj) 



C A a s (iJ,) 



\n(l-r 2 )S'(£). (70) 



Note that Cf(£,r,n) must have mass dimension —2. At 
0(a s ) the simplest potential solution has the dependence 
ln(/x//t) $'(£), but by dimensional analysis the only pos- 
sibility for k is £ which leads to a singular result. The 
correct solution is 



S'(£)- 



C A a s (ji) 2 d 
In 1-r — 

5'{£)d\{r)+0{a 2 s ), 











.£- 


♦) 



(71) 
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which can be deduced since the derivative of the plus 
function has the right dimension and has the required 
logarithmic scale dependence 



_d_ d_l r 

d/i d£ fj, 



S'(£) 



(72) 



In this way, the plus function term in the Wilson coeffi- 
cient exactly compensates for the first order a s (n) ln(/i) 
dependence in f2i(r, fi). Note that 



d_l 



7'~ 




1 


\K 2 ] 


.£. 


+ 


~V 2 


.£ 2 . 



1 

/' 



6(e) 



(73) 



where the H — \~ distribution induces two subtractions 
about £ = and is defined so that its zeroth and first 
moments integrate to zero for the limits £j [i € [0, 1]. 

The function d\ (r) in Eq. (jTTj) is also a perturbatively 
computable contribution to the Wilson coefficient. The 
matching calculation for this term involves considering 
the difference between renormalizcd Feynman diagrams 
for the full theory soft function matrix element 



S e {l) = (Q\Y\Y^{£- 



Qe)Y n Yn\0) 



(74) 



and for the low-energy matrix elements describing 
S!i(r, fi). A complete one- loop calculation of d\(r) is be- 
yond the scope of our work. In App.[D]we carry out this 
matching procedure for thrust in order to directly derive 
the term that involves the derivative of the plus-function 
shown in Eq. (|71|) . Many of the complications required 
to derive d\ (r) do not enter for this term. 

Next consider the impact of Cf (£, r, fi) on the distribu- 
tion and first moment event shape observables discussed 
in Sec. Ill CI in order to determine the appropriate scale 
fi where large logs are minimized. For the differential 
distribution we find 



da da 1 

de de C 



da ( 

d 

x c e g e (r) ttx(r, 



Q 



Q 

:Q d£ 



d*_a 
de 



d 2 a 



, (75) 



where the nonperturbativc parameter Of' n (/x) is given 
in Eq. dHZJ) , and 



^i I (/ i )= /drdf(r)ce5 e (r)Oi(r, /i). 



(76) 



The explicit ln^/Qe) in Eq. (|75|) implies that for the 
distribution, the appropriate scale to run the power cor- 
rection to is \x = Qe. 



For the first moment we find 



ul d l( e -^)F e (£) 



a de 



Q 



(77) 



= de d£ 



£lf(li) a s (/j,) n e { d nl'^ifi) C A a s (fi) 



i pert 



T 1 da 
x / de- — (e) 

a de 



In 



7T Q 
^t?(^max 



Q 



c(^max ^) 



where the notation (e) rt is defined in Eq. (|20|) . Since 
the perturbative moments generate a rapidly conver- 
gent series, we can expand the perturbative coefficient 



of f2^' ln (/x) in (e/e m ) to obtain 



a del^(e) 
a de 



= In 



In 





e 2 


max ^0 




( e )pcrt 


( e2 ) ♦ 

\ / pert ! 


^max 


LC, max 



(78) 



Thus for the first moment the appropriate scale to run 
the power correction to is ji — Qe max . 



C. Orthogonal Basis for f2i(r, /x) 



In Sec. IIVD1 we showed that the power corrections for 
observables in different universality classes could be ap- 
proximately related by expanding out g e (r) and ^i(r, fi) 
in a suitable basis. Since the appropriate scale (i for fii 
is Q dependent, this is true if all measurements are per- 
formed at a single Q. From Sees. IV Al and fVBl we know 
that the leading power correction has nontrivial Q depen- 
dence, and we would like to incorporate this information 
in our description of ili(r). 

Immediately from Eq. we see that fli(r, /i) will di- 
verge as r — > 1. Even if f2i(r, /xo) is regular at some scale 
/io, it will quickly develop a singularity at r = 1 as fx 
evolves. This singularity can be physically interpreted as 
the propensity of Wilson lines to emit soft masslcss parti- 
cles. Of course, this singularity is still square-integrablc, 
and thus the power correction is still well-dcfincd03 How- 
ever the singularity at r = 1 means that if the Legen- 
dre polynomial basis in Sec. IIVDI is used at a scale fio 
to define a basis for power corrections fij (/^o), then it 
converges very slowly when trying to use these same pa- 
rameters to describe power corrections for scales fi fiQ 
or fi <C Ho- 



12 For evolution over a large enough Q range, the singularity may 
no longer be integrable, and may turn into a distribution, but 
we do not encounter this subtlety in our analysis. 
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Looking at the a s expansion of the running formula 
in Eq. ([6"B")). we see that to describe the power correction 
for a range of Q values, we must find a suitable basis to 
describe not only g e (r) and fli(r, /io) but also g e (y) ln(l — 
r 2 ). Since ln(l — r 2 ) is unbounded at r = 1 (but still 
square integrable), we will use the freedom to introduce 
additional square-integrable basis elements (1 — r)~ k for 
suitable values < k < 1/2. As long as we are content 
to work with a finite number of basis elements, we can 
always make such a basis orthonormal via the Gram- 
Schmidt procedure. Adding these additional functional 
forms will in general yield an over-complete basis, but 
this is not an issue in practice since we will only ever 
consider a finite number of basis elements. 

The situation becomes a bit more complicated if we 
consider the full running in Eq. (|64[) . In general, any 
finite basis we choose to describe S7i (r, /j.) at one value 
of Q will not provide a good description at a different 
value of Q due to the running. In particular, a basis that 
is orthonormal at one value of /io will no longer be or- 
thonormal at another scale. Instead of trying to find a 
basis that works for any Q, we instead choose a scale [io 
at which we model 51i(r, /to), and then evolve according 
to Eq. (|64|) . We then fit for the basis coefficients by us- 
ing information at different values of Q. This procedure 
is philosophically the same as the procedure used to de- 
termine parton distribution functions (PDFs) , where the 
PDFs are modeled at a low scale and then evolved to 
higher Q values via the DGLAP equations [75l - l77j . 

D. Comparison to Monte Carlo 

We now show that power correction universality and 
running is exhibited by two widely used Monte Carlo pro- 
grams: Pythia 8.162 [78| and Herwig++ 2.6.0 [79|. The 
hadronization model in Pythia 8 is based on string frag- 
mentation while Hcrwig++ is based on cluster fragmen- 
tation!^ The default hadronization parameters in both 
programs have been tuned to reproduce LEP e + e~ event 
shapes at the Z pole. We will consider e + e~ — > hadrons 
at various Q values, turning off initial-state electromag- 
netic radiation to avoid the radiative return process. 

Our study will use the generalized angulari ties Tf n a) 
defined in Eq. (|37|) . Via the arguments in Sec. HVBl we 
know that the power corrections for different values of a 
are related via the c a = 2/(1 — a) coefficients in Eq. (f39|) . 

n T 1 { - ) (p)=c a W(p), (79) 

where f2™ (/•*) i s the universal power correction for the r n 
class given in Eq. (|42l) . 



The two programs also have different showering models, with 
Pythia 8 using a pj_ -ordered shower and HerwigH — h using an 
angular-ordered shower. Since the showers are evolved down to 
the non-perturbative scale, part of the renormalization group 
evolution of the power correction may be captured by the show- 
ering algorithm, and not just the hadronization model. 



From Eq. (|2"Tj) . we know that the leading power cor- 
rection shifts the moments of event shapes, so we can 
extract the universal power corrections for the general- 
ized angularities via the first moment 

(Wj) = T( Q < T ('^>) - Q ( T ('^))port) . ( 80 ) 

up to small higher-order corrections. Since the maxi- 
mum value of thrust (a — 0) is 1/2, we take [aq = Q/2 to 
avoid having large logs in these higher order corrections, 
which were displayed above in Eq. ([75| F^I The pertur- 
bative moment ( r (n.a)) pcrt lS the same for event shapes 
with a common r — > 1 limit. To form a combination that 
is sensitive to the power corrections in the Monte Carlo 
programs without having to know about their perturba- 
tive contributions, we consider the difference T(o. a )~ r (n.a) 
which compares the same value of a at two different val- 
ues of n. For this combination we have 

Q?( MQ ) - fi?( Mo ) = 9- ( (r C0 , o) ) - (r (n>a) ». (81) 

Ca 

Note that this difference is also independent of the ad- 
ditive scheme change that removes the renormalon from 
fif(/x), given by Eq. (|55|) with Eq. (f58|) . Therefore we 
emphasize that our analysis in this section only probes 
the running in Eq. ([B"3")l and not the devolution [5^, [6(3] 
associated with ill(R,fi). 

In Fig. Uwe plot fi? - fi? for a = 0, -1, -2, n = 1, 2, 3, 
and Q values ranging from 20 GeV to 200 TeV. At a 
fixed scale Q, the power corrections are independent of a 
with at most 5% variations, thus demonstrating the an- 
ticipated universalitvP'l Note that both programs were 
tuned to LEP Z pole and low energy data, so it is not 
surprising that they have the same power corrections at 
Q ~ mz- More interestingly, both programs show loga- 
rithmic growth in Q for the power correction, as expected 
from our results in Sec. IV Al Numerically, this growth is 
consistent with the form (In Q) A /Q found in Ref. [l|, and 
the exponent A ~ 4Ca/ /3q ~ 1.5 is presumably related to 
the exponent in Eq. (|64[) . A more concrete comparison is 
difficult since the analysis in Ref. [l[ effectively expands 
about r = 1, and parametrizes the extra resulting loga- 
rithmic singularity by a 1u(/x/Aqcd) factor that cancels 
an a s ((i). 

To show the importance of resummation, we fit for 
the functional form of f2(r, fio). The solid lines in Fig. [5] 
have the full running in Eq. (|64[) . while the dashed lines 
correspond to the expansion in Eq. (|6"6"|) . These curves 
were obtained following the procedure of Sec. lVC| where 
f2(r, no) is modeled using the three basis functions 

{l,r,(l-r)- 1/A }, (82) 



14 The angularities with a 5^ 1 have different values for e max , but 
this is an 0(1) change to the fi scale, and hence not relevant. 

15 The leading violation of universality can be attributed to differ- 
ent matching coefficients d\(r) in Eq. (1716 - 
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FIG. 8. Universal power corrections extracted from 
Pythia 8.162 and HerwigH — h 2.6.0 for the generalized angular- 
ities T( n ,a) ■ We use the measured power corrections in the first 
two plots to fit for 0(r, (Ms) at Q = 100 GeV, using Eq. ([511 to 
evolve to different Q scales. The curves in third plot are then 
predicted. The dashed curves show the approximate formula 
in Eq. <[EB]l. 



suitable orthonormalized. The inclusion of (1 — 7-) -1 / 4 
is needed to capture the (intcgrablc) peak of Q(r, fio) at 
r = 1, though other choices (1 — r) k give comparable 
results. We apply the fit form at Q = 100 GeV and use 
Eq. to determine f2(r, /i) over the whole Q rangeP^I 
To show that this framework has some predictive power, 





Q = 100 GeV 
Pythia 8 Herwig++ 


Q = 10 4 GeV 
Pythia 8 Herwig++ 


ol - o\ 


0.7 GeV 0.6 GeV 
0.9 GeV 0.5 GeV 


1.3 GeV 0.9 GeV 

2.4 GeV 1.1 GeV 


ol - nl 

O"' - 0{ 


1.1 GeV 1.0 GeV 
1.8 GeV 1.0 GeV 


2.3 GeV 1.6 GeV 
4.6 GeV 2.0 GeV 


ol - ol 


1.5 GeV 1.3 GeV 

2.6 GeV 1.3 GeV 


3.2 GeV 2.1 GeV 
6.6 GeV 2.8 GeV 



TABLE VII. Power correction differences extracted from the 
fits in Fig. [8] The slope parameter O"' ln is defined in Eq. (|83|) . 
These values have 10% to 20% uncertainties from the choice 
of functional fit form. 



we fit f2(r, fig) using information from n = 1 and n = 2 
over the whole Q range, and then extrapolate to n = 3. 

Because we are measuring the difference T(o )0 ) — T (n,a)j 
we are relatively insensitive to the functional form of 
n(r, no) near r = 1. We are effectively measuring event 
shapes in the universality class g e (r) = 1 — r n where 
p e (l) = 0, and therefore our extraction of the raw power 
correction f2" has large uncertainties. On the other hand, 
the extraction of power correction differences is stable at 
the 10% to 20% level as the fit form is adjusted, and these 
differences arc shown in Table IVIII for the basis choice in 
Eq. (|82]). 

We also show the best fit values for f2?' ln , defined as 



^i' ln Oo) 



dr ln(l-r 2 )r"Qi(r,Mo), (83) 



which is a hadronic parameter related to the slope of the 
expanded running in Eq. (|66|) . We see that Pythia 8 
and HerwigH — h have similar power corrections f2™ at 
Q = 100 GeV, but the slopes f2™' 11 are larger in Pythia 8, 
leading to larger values of f2™ at Q = 10 4 GeV. This slope 
parameter is interesting since it provides an example of 
a hadronization effect that can only be accurately deter- 
mined using data at multiple Q values. Parameters of 
this type presumably dominate the uncertainty one has 
when describing hadronization effects in high energy data 
using Monte Carlo models that were only tuned at much 
lower energies. 

Finally, we remark that one strategy to extract f2" di- 
rectly from these Monte Carlo programs would be to turn 
off the hadronization model and calculate ( T (ri,a)) pcrt 
from the parton shower alone. However, there is no guar- 



As mentioned below Eq. H64I) , a more natural strategy would be 
to apply the fit form at fio = 2 GeV. Because of large range 
of scales in Fig. [8] resummation is always important with that 



choice. By using /^o = (100 GeV)/2 we are using the same scale 
where Monte Carlos have been tuned, and we can also better 
highlight the difference between the full and expanded running. 
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antee of a one-to-one map between hadronization model- 
ing and operator-derived power corrections. In the con- 
text of a Monte Carlo program, the perturbative par- 
ton shower is first evolved to the shower cutoff of order 
1 GeV > Aqcd before applying the hadronization model, 
whereas from Sec. IV Bl the natural scale to evaluate the 
power correction is Qe (for the distribution) or Qe max (for 
the first moment). For this reason, there is an ambigu- 
ous separation between perturbative parton shower evo- 
lution and nonperturbative hadronization modeling, and 
there is no guarantee that hadronization models by them- 
selves will respect the renormalization group evolution of 
Eq. flUEZI 

It would be interesting to understand to what 
extent parton shower evolution can mimic Eq. (|64]h and 
whether there are ways to adjust hadronization models 
to satisfy the renormalization group properties expected 
of power corrections. Ultimately, one would want to test 
the power correction evolution by performing event shape 
measurements at high Q. 



VI. CONCLUSIONS 

In this paper, we revisited the important issue of power 
corrections for e + e~ dijet event shapes. By casting the 
leading power correction in terms of matrix elements of a 
transverse velocity operator, we are able to robustly treat 
the effect of hadron masses. Depending on the measure- 
ment scheme, event shapes fall into different universal- 
ity classes that share a universal power correction Of" . 
Moreover, these nonperturbative matrix elements have 
perturbatively calculable anomalous dimensions, which 
introduce additional dependence on the scale Q of the 
hard collision. 

Since Monte Carlo programs play such a key role in 
LHC data analysis, it is satisfying to see that both univer- 
sality and Q-evolution are exhibited by the hadronization 
models of Pythia 8 and Hcrwig++, albeit with different 
choices of the nonperturbative matrix elements. An in- 
teresting difference is in their values for hadronic slope 
parameters that play an important role in the extrapola- 
tion to high energies of hadronization effects that are fit 
at low energies. 

Our study motivates a reanalysis of e + e~ event shape 
data with a more explicit treatment of hadron mass ef- 
fects. As an exercise, we have studied the effect of includ- 
ing the anomalous dimension of the leading power cor- 
rection on the determination of a s (mz) from the thrust 
distribution. We have repeated the analysis of Ref. [3(3] 
at N 3 LL + O(al), using the same data set and procedure, 
but replacing the power correction to include a logarith- 



17 For example, if we were to extract fij from Pythia 8 by turning 
hadronization on and off, we would find almost no running with 
Q (i.e. Q°' ln ~ 0). 



mic slope term 

2 0l _> 2 0l _ 2^}£a nj- log , (84 ) 

where (J- s (t) ~ Qt. With the current experimental data, 
Oi and Si^ are highly correlated and cannot be deter- 
mined simultaneously. Therefore we plug in the estimate 
SI 1 / 1 = ±Oi, and repeat the fit for a s (mz) and fii. We 
find that the effect of this term on a s (mz) is ±0.0005, 
which is roughly half of the total uncertainty ±0.0011 
found in Ref. [30|. For ft x the effect is T 0.03 GeV, which 
is comparable to the previous uncertainty of ± 0.05 GeV. 

If LEP data is successfully preserved [8(1 HH, then 
one could compare different mass scheme to better sepa- 
rate perturbative physics from nonperturbative physics. 
Such studies would be interesting in their own right, but 
would also give additional input for tuning Monte Carlo 
hadronization models. 

An obvious generalization of this study is to go beyond 
dijet event shapes and consider shape functions for more 
than two Wilson lines. This is relevant not only for mul- 
tijet studies at e + e~ colliders but also for treating the 
beam directions in hadron colliders like the LHC. For ex- 
ample, the event shape iV-jettiness [(34j| is a convenient 
variable to define exclusive AT-jet cross sections at the 
LHC, and its shape function involves 2 + N Wilson lines. 
The anomalous dimension calculation for multiple Wil- 
son lines is technically more challenging but conceptually 
similar to the calculation presented here. In particular, 
the same £T{r,y,t) can be considered for this analysis 
(where here we add i in order to emphasize that rapidity 
y is defined with respect to the axis t). On the other 
hand, it is not clear which aspects of universality will 
carry over to the multijet case, since universality of dijet 
power corrections relied crucially on longitudinal boost 
invariance. 

Finally, our study sheds light on recent studies of jet 
substructure at the LHC. The jet shape A^-subjettiness 
was introduced in Refs. [H, HH (see also Ref. [Hj]) as a 
complement to the event shape iV-jettiness. The ratio 
of 2-subjettiness to 1-subjettiness (t 2 /i) can be used to 
distinguish boosted W/Z bosons from the background of 
ordinary quark and gluon jets. From Table [J we see that 
2-jettiness and thrust are closely related, and Ref. (85[ 
used this fact to perform a precision calculation of t 2 /\ for 
boosted W/Z bosons at the LHC by recycling the known 
resummation for thrust in e + e~ — > hadrons. However, 
Table IVl shows that the original definitions of 2-jettiness 
and thrust are not in the same universality class, which 
explains why Ref. [85| required a value of the power cor- 
rection <f>i = ill that was roughly a factor of two bigger 
than the value for Sl\ obtained in Refs. [2^, [3(J. To a 
good approximation, the power correction for boosted 
W/Z bosons generates a shift of the t 2 /i distribution by 
7T times SI] [85[ , making it all the more crucial to choose 
the proper power correction. Thus a proper treatment 
of hadron masses and power corrections will be essential 
for precision jet substructure studies at the LHC. 
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Appendix A: Derivation of the Transverse Velocity 
Operator 

In this appendix we show that the transverse veloc- 
ity operator £T(r,y) can be expressed in terms of the 
energy- momentum tensor as in Eq. (|26[) . Our analysis is 
analogous to that for £t(v) m R- e f- H3j where it was per- 
formed for scalar and spin- 1/2 hadrons, as well as that 
of Ref. [61| . We will carry out our proof for scalar fields. 

Consider the energy-momentum tensor of a free scalar 
particle with mass m (we will see below that interaction 
terms are suppressed): 



(Al) 



where the plane wave expansion for the scalar field is 



where E p = \/ p 2 + m 2 . We will use the stationary phase 
approximation 



lim dxf(x)e lk9{x) 



(A3) 



2n 



k\g"(x )\ 



f{xo) e 



ikg(x ) e i f sign[g"Oo)] 



where g'(xo) = 0. For this formula to be applicable, g(x) 
must attain a minimum or a maximum in the range of 
integration. 

After plugging the plane wave expansion of the free 
field 4> in Eq. (|A2[) into the energy-momentum tensor in 
Eq. (|A1|) . one can perform all the angular integrations 
using Eq. (|A3|) . obtaining 



lim R 3 h i T 0i {R 1 Rvh) 

R— ¥00 



a iR(vp-E p +vq-E q ) 



v d f dpdqpq 2 f 

= lim H / — ; 7-r— — - (Iff. dff- 

R^oc J A(2TT) 4 E q V 2 V P " 9 " 

+ a P* 4* e iR(vp - E "- vq+E * } + (h.c. and p , 

where p„ = ph and qa = qh. The p and q integrals can 
be performed, again using Eq. (|A3|) . to yield 



lim R hi T (R,v Rh) 

R— >oc 



m v 



4(2^) 3 (l-^a); 



■ (a p at + at a^j , (A5) 



where p 



h. Note that terms involving two cre- 



ation or two annihilation operators drop out, since they 
vanish as 1/R when integrated against any function of v 
(in particular our g e (r) in Eq. (|36p V At this point we can 
also see why interaction terms in the energy-momentum 
tensor can be neglected. Such terms involve additional 
fields and therefore additional integrations when they are 
expanded in plane waves. These additional integrals can 
be performed using the stationary phase approximation 
and vanish as R — > oo due to the presence of additional 
powers of 1 /R. 

After normal ordering, Eq. (|A5[) can be written as 



lim R h.,, T (R, v Rh) 

R^oo 



(A6) 



d 3 P t E p f \p\ 
(27r) 3 2£ p 



4*£'(-|!)<"tf-*>- 



Using the fact that E p = (p 1 - cosh rj) jv and 

d<pS 2 (p — h) = cosh 2 ^ 8(i] — 7] p ) . 

o 



(A7) 



we obtain the operator 



&r(v,rj) 



lim R 



2n 

3 / J ± * rrtOi 



cosh T] 



hiT Ul (R,Rvh), (A8) 



which is differential in velocity and pseudo-rapidity and 
satisfies: 

e T (v,r,)\X) = J2pl5(v-v i )5(r,-r H )\X) . (A9) 

Note that if we integrate this operator over < v < 1, 
we recover the expression in Ref. [6l| for the energy flow 
operator £t{v)- 

Finally, to obtain from £t{v,t]) the desired £T(r,y) 
that satisfies Eq. (f25j), one needs to multiply by the Ja- 
cobian factor 



d(v,rj) sech 2 y 
d(r,y) r 



(A10) 



and include a factor 1 jr to convert p to m , yielding 



sech y 



£r(v(r,y),T}(r,y)) 



(All) 



(A4) 

which agrees with Eq. (|26l) . 
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Appendix B: Renormalon Computation for Generic 
Event Shape 

Here we show that the definition 5 e {R, fi) = 
(c e /c e i)S e i (R, fi) in Eq. ([55} yields a perturbative cross 
section (T e (x) in Eq. (|56|) that is independent of the 
leading Aqcd renormalon when probed by a standard 
fcrmion bubble chain. Since the renormalon cancels be- 
tween the MS series a e {x) and r2f(/i), this implies that 
£ll(R, n) is also free of the Aqcd renormalon. 

The Aqcd renormalon corresponds to a u = 1/2 pole 
in the Borcl transform. For a function f(a s ) that is an 
infinite series in a s {fi), the Borel transform B[f]{u) is 
obtained by replacing 



/ /3 q^h " +1 ^ t£ 
I 4i / n! 



(Bl) 



Following Rcf. [57| we make use of the fact that the per- 
turbative soft function carries the leading renormalon, 
and hence carry out our computation for Sf cvt (x, /i) 
rather than the cross section <r(;r)0£] Since the soft func- 
tion obeys non-Abelian exponentiation [8(1 [87} it is useful 
to write the perturbative scheme change in Eq. (|5o) as 



lnSP elt (x,^) =lnSP crt {x,Li)-ix$e(R,v), (B2) 

and then demonstrate that In S% eTt (x, //) does not have a 
u = 1/2 pole. 

The use of the soft function allows us to perform the 
bubble chain analysis for an arbitrary event shape spec- 
ified by f e (r,y) and Eq. (|9]). To study the first contribu- 
tion to the u = 1/2 pole, we can work in d = 4 dimensions 
and we only need to dress a single real gluon with a bub- 
ble chain. We parametrize the gluon phase space with 
p ± and y, 



d 3 p 



_ dy d 2 pj_ 
{2ir) 3 2E p ~ 4tt (2tt) 2 



(B3) 



Since the final state gluon is on-shell, we have r = 1. For 
the event shape e, the Fourier transform gives 



6 ~ Q p± ^ 1,y ^ 



(B4) 



Taking the sum of all dressed real radiation diagrams 
with a single gluon and swapping rif — > — 3/3o/2, we find 
the Borel transform 



B[ln5 bubblcs ( ^ )](u) 

8C F ( M 2 e 5/ T 



A,r(i+u)r(i-«) J_ 

8C F (M 2 e 5 / 3 ) M 

/3 r(i+u)r(i-u) 



/-\-oo poo 
dy / dp ± p~ 
-oo JO 
hoo 

dZ// e (l,2/) 



(B5) 

-ixp±f e {X,y) 



-\-oo poo 

2u I dh h- 1 - 2u "- ixh 



Note that x is a dimensionlcss variable in a-(x) but is a variable 
with mass dimension —1 in Sp CTt (x, fi). 



^ r(i+u)r(i-«) 



Here /3o = 11CU/3— 2n//3, and in the second equality we 
used the change of variables h = p±f e (l,y). Expanding 
about u = 1/2 and using J dy / e (l, y) = c e , we arrive at 
the final expression for the u = 1/2 pole 



B[ln^ ubblcs (x ;At )] (u) = Ce 



8CVe 5 / 6 



7T/3 (u - a) 



3- (ia?/i) . (B6) 



Here (ix) corresponds to the S'(£) present in Eq. (|T6|. 

Using Eq. (|B6|) we can compute the Borel transform 
of the subtraction series 5 e i (R, fi) for the reference event 
shape e', which is defined by Eq. (|58[) . We find 



B[S e ,{R,[i)](u) = c e 



8C F e 5 / 6 



/' ■ 



(B7) 



Finally, computing the leading renormalon ambiguity in 
ln»Sp 0rt , using Eq. (|55|) to define 8 e (R,n), we find 

B[mSr rt (^)]H (B8) 
= B [ In 5 e port (x, /i)] (it) — ix B [6 e (R,n)](u) 

= B [ In 5P° rt (x, //)] («) — ix — B [5 e , (R, /*)] («) 

o 



as promised. Note the importance of using the same 
scale for the perturbative soft function 5'P ert (x, fi) and 
its subtractions S e {R, jti). 

As a final comment, we remark that the renormalon 
analysis in this appendix takes r = 1 and hence does 
not fully probe infrared effects that depend on hadron 



Appendix C: One-Loop Anomalous Dimension 

In this appendix we provide details on the calcula- 
tion which yields the anomalous dimension formula in 
Eq. (|61[) . The integrals involved in determining the one- 
loop anomalous dimension of f2i(r, fi) from Figs. 2]-[7]are 
somewhat different from the phase space integrals for 
QCD gluons attached to eikonal lines, which occur for 
the leading power perturbative soft function calculation. 
In particular the amplitudes are similar to those occur- 
ring in recent two- loop soft function calculations (88l490| , 
but a different measurement is made. 

As explained in Sec. [V] one of the cut gluons corre- 
sponds to a massive adjoint source field, and £t{t, y) acts 
on this object. We will call the momentum of this source 
q^, where q 2 = m 2 ^ 0. The remaining gluon lines are 
standard masslcss QCD gluons, and we will denote the 
momentum of virtual loop integrals by k, and the mo- 
menta of real gluon radiation by p, where p 2 = 0. 
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The phase space integral over q for the source is com- 
pletely fixed by the measurement, up to one trivial an- 
gular integral for <\> q in the transverse plane. Taking the 
three phase space variables to be q + , q~ , and <f> q the for- 
mer two are fixed by the 8- functions from Eq. (|25|) : 



mjd(r - r q )S(y - y q ) 
2 m 2 r 



(CI) 



m 



<i (1 _ r 2)2 



51 q + - 



me" 



vT 



me 



vT 



where mr- 



q±/rrig, and y 9 



1/2 \i\(q + / q~). For notational convenience wc define an 
object $(r, y) that contains common factors associated 
with the source that appear in all Feynman diagrams, 



$(r, y) 



l6Tra s C F f d n q m£ 6{r-r g ) 5{y-y q ) 
(2tt)« J 2E q 



q+q- 



(C2) 



where n = d — 1 = 3 — 2e. Here $(r, y) should be consid- 
ered to be an operator that can act on additional q + and 
q~ dependence in loop and phase space integrals, and 
which replaces q ± by the functions of m, r, y occurring 
in the 5- functions in Eq. (|C1[) . <&(r, y) is normalized so 
that acting on unity with n = 3 the integral in Eq. (|C2j) 
yields the tree level result in Eq. ([50)) . 

The general strategy to compute the anomalous di- 
mension is to reduce each graph to a set of master inte- 
grals. We will always partial fraction eikonal propagators 
and shift numerators to obtain a set of integrals that in- 
volve only one p + (k + ) and/or one p~ (k~) in a denom- 
inator. To regulate potential IR singularities we shift 
the eikonal propagators by taking ^ — > + A ni n or 
p ± — > p ± + A n ,n. We will treat the A„ jfi as infinitesimal 
IR regulators, which are expanded whenever possible. 

We start with graphs which involve independent emis- 
sion of gluons in Figs. [5] and El A useful identity between 
real emission phase space and virtual integrals is 



h(A,B) =/2 



2\p\{2ir) d - 1 (p++A) (p~+B) 
d d k 1 11 



(27r) d (fc+ + A) (k~ + B) k 2 + iO 

where p, = /ie 7E /(47r)03 Using Eq. ([C3| it is easy to 
determine that the purely Abelian terms proportional to 
Cp vanish (virtual graphs cancel real radiation graphs). 
In addition in the sum of non- Abelian contributions from 
the graphs in Fig. El there are no IR divergences regulated 



19 Eq. l)C3[l has IR divergences regulated by e for either A or B 
zero (strictly speaking UV divergences cancel IR divergences in 
pure dimensional rcgularization, and the integral is zero). The 
equality between real and virtual integrals is still valid for these 
cases. 



by e or A„ :S , so all 1/e's correspond to UV divergences. 
The result for these non-Abelian independent emission 
contributions is 



EFT 



InasCA^ir, y)h(q + ,q~ 



(C4) 



2tt 



<P{r,y) 



r 1 



In 



<T<f 



7i~ 

T 



+ 



Note that we refer to the results of this section as effective 
field theory (EFT) contributions since we are performing 
calculations in a theory where the soft perturbative scale 
(MS ~ Qe) has been integrated out. 

Graphs involving the triple gluon vertex are more in- 
volved. Here individual virtual radiation graphs have 
an imaginary part, which arises from the fact that 
q 2 = m 2 > while all virtual particles have massless 
propagators. Nevertheless the virtual diagrams can al- 
ways be paired with a complex conjugate so one only 
needs the real parts. To determine the anomalous dimen- 
sion for Qi(r,[i), we only need the UV-divergent terms. 
The finite terms would be necessary for the matching 
computation of d\{r) in Eq. (|7ip . which is not our goal 
here. We will therefore focus on graphs which contain 
1/e UV divergences and \i dependence. The ^-dependent 
terms will be needed for App. [D] 

It is straightforward to verify that the sum of double 
cut vacuum polarization graphs in Fig. U does not have a 
1/e UV divergence. They also do not require e to regulate 
IR divergences, and hence have no explicit fi dependence. 
Thus they do not contribute to our calculation here. 

The remaining triple gluon vertex diagrams shown in 
Fig. [7] involve single cut virtual graphs and double cut 
real emission graphs. The only required UV divergent 
and /i-dependent loop integral is 



I 2 {A 7 q+ 7 m) 
= Re 



(C5) 



ifi 2 ' 



d d k 



(2n) d (k+ + A)k 2 + iO (k - q) 2 + iO 



(4 

The result is IR safe for all the cases we need (A ^ and 
A ^ -q + ), yielding 



h{A, q + , m) = - 



rl 



(4tt)V 



(C6) 



where the + . . . refers to UV finite and /^-independent 
terms. The real radiation master integral is 



h{A,q + ,m) =/i 2 



id-l 



1 



1 



2\p |(27r) d - 1 (p+ + A) {2p ■ q + m 2 ) 



(47T) 
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which is IR safe for all cases we need (A 0), yielding 
1 1 f 1 1 



h(A,q + ,m) 



(4tt) 2 q+ 2e 2 e 



In 



(C7) 



1 ffn 
- - 2 In — 
e V 



Here all 1/e poles are UV divergences, and again the 
+ . . . represent UV finite and jU-indcpcndent terms. We 
can always group terms with IR divergences as A„ ifi — > 
into the IR finite combination 

I 23 (q+,m) = h(A,q + ,m) + I 2 (A,q + ,m). (C8) 

The result is independent of A 



/ 23 (? + ,m 2 ) 



1 1 

(4tt) 2 ^+ 



1 

2c 2 



1 / 111 
— In I — 

e V /i 
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(C9) 



EFT 



The sum of all diagrams in Fig. [7] is 

Ana s CA^{r, y) q~ h?,{q~ , m) + q + hz{q + ,m) 
+ q~h{q~,q~,m) + q + I 3 {q + ,q + ,m) + ... 

$(r,y) 



C A a s 



2tt 



(CIO) 



Finally, adding all one- loop diagrams, the 1/e 2 terms 
vanish, and we get the following final result for the UV- 
divergent and /^-dependent terms from the graphs in 

Fig.m 

+ .... (Cll) 

Due to $(r, y), we have \n(q + q~ /m 2 ) = — ln(l — r 2 ), and 
this result reproduces the 1/e term shown in Eq. (|0 1 [) I 2 "! 



20 This same result can also be obtained setting A n> ^ = from the 
start and using dimensional regularization for the IR divergences. 
In this case we would use 

h(0,q + ,m 2 ) = -I 2 (-q + ,q + ,m 2 ) 

_i r(e)r(-e)r(i - e) 1 



and 



l 3 {0,q+,m) = -I 3 (q+,q+,m) 



(4tt) 2 T(l - 2e) 

r(2e)r(-e) 1 



cos(e7r) ( e lE 
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-1 1 

(4tt) 2 ~q+ 
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(47r) 2 q+ V m 2 / 

-ln(-) +ln 2 (-) + 

c V fj,J \ 24 



In this case the sum of all diagrams is again IR finite and again 
reproduces Eq. <C1U . 



The expression in Eq. (|C11|) corresponds to a bare re- 
sult in the EFT. The renormalized result Mf^J p is simply 
obtained in minimal subtraction by canceling the UV- 
singularity with a countcrterm Z(r, fi, e), and will be used 
in App.[D] For that purpose, we also note that the sum of 
independent emission and triple-gluon virtual graphs is 
UV finite, so the entire contribution shown in Eq. (|Cllj) 
comes from the real emission diagrams. 



Appendix D: One-Loop Matching for Thrust 

In this appendix we consider the one-loop matching 
computation that determines Cf (£, r, fi) in Eq. for 
the case where e is thrust. Only ^i-dependent terms will 
be considered since our goal is to see how the function 
d/d£[l/fi(/j,/£) + ] arises from the matching computation. 

A matching computation is performed by considering 
the difference of renormalized full theory and renormal- 
ized effective theory matrix elements, which are calcu- 
lated with precisely the same infrared regulator(s). For 
our computation the full theory corresponds to matrix 
elements of the soft function. For an event shape e, it is 

iTTt 



S e (£) = ( I YlY^S(e - Qe)Y n Y n \ ) (Dl) 

= I YIy 7 ] \X)8[£-Qe(X)]{X\ Y n Y n | ) , 
x 

where the sum is over all possible intermediate states 
and includes also integrals over phase space. The EFT 
for this computation corresponds to the field theory ob- 
tained by integrating out the scale £, and involves ma- 
trix elements like £li(r) in Eq. (|32|). In this language the 
anomalous dimension calculation in App. [Cl corresponded 
to finding the UV countcrterm for the fli(r) matrix ele- 
ment in the EFT. For a matching computation that only 
considers /x-dependent terms, the required renormalized 
EFT matrix element corresponds to the ln(m///) term in 
Eq. (jClip . This result is independent of the event shape 
e, in contrast to the d\{r) term in Eq. (|71[) which is event 
shape specific. 

Here we perform a computation of the corresponding 
full theory matrix elements in Eq. (jDll) with the same 
IR regulators, which includes the use of the source J^ A 



field of momentum q^ where q 



Since some parts 



of this computation depend on the choice of e, we will 
restrict ourselves to a computation for thrust. We will 
use the same notation as App.[C]ior loop and phase space 
integrals. Unlike the EFT, the full theory results involve 
hierarchical scales, £ 3> q^ ~ m, and hence the final full 
theory result must be expanded before subtracting the 
EFT result. 

At the order of our calculation, we can split the 
in Eq. (|D1[) into terms with no source term (which gives 
rise to the purely perturbative soft function), and terms 
with one source term. In general the expansion £ q 
should only be performed after carrying out the full the- 
ory loop and phase space integrals, but in cases where 
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the expansion and integration commute we can do them 
in either order. One example where this is useful is in the 
measurement function 8(1 — Qe(X)). Denoting by e(q) 
the contribution from the source, and e(j>i) the contribu- 
tion from all other real radiation gluons we can expand 
I r~j Qe(pi) 3> Qe(q), and keeping the first two terms 

6[t-Qe( Pi )-Qe(q)} (D2) 
= 5[£ - Qeipi)] - Qe(q) S'[£ - Qe( Pt )] . 

For the corresponding terms in S e (£) this yields 



{pi},q 



x (Q\YlY 7 l\p l ,q)(p. l ,q\Y n Y n \0). 

The first term corresponds to the leading power pcrtur- 
bative soft function and the second term provides the full 
theory contribution to the matching we are interested in. 
The analog of Eq. (|D2[) for the EFT computation is 



6[£ - Qeipi) - Qe(q)} = S(£) - 5'{£) [Qe(q) 



Qe(pi)] , 
(D4) 



where the term Qe(pi) is scaleless and vanishes. From 
this result we see that the EFT contribution to the 
matching is proportional to —d'(£). 

It is easy to check that the matching is simple for an 
Abelian theory. Due to the exponentiation and factoriza- 
tion properties of the Abelian cikonal matrix elements, 
we obtain 

- ]T (0 1 Ft y B t | Pl , q ) 6' [£ - Qe(pi)} Qe{q) ( Pi ,q\ Y n Y n |0> 
= -^r rt W(0|Fiy„t(Qe)F„F« |0) 

= ~rt^. (D5) 

This result holds even if we consider including more 
than one source term. In Ref. [3(| it was assumed 
that Eq. (|D5|) also encoded all non-Abelian contribu- 
tions. While these non-Abelian contributions are indeed 
present, Eq. (|7ip implies that in general there are addi- 
tional non-Abelian corrections from the +- function and 
df(r) terms. The flaw in the argument in Appendix B 
of Ref. |3(j is that the dimension- 1 operator £t(v) 1S n °t 
unique, since there exists an entire family of operators 
£T{y,r) parametrized by r. 

To carry out the full non-Abelian calculation for the 
second term in Eq. (|D3|) . we use Eqs. (|25|) and (|27|) to 
decompose Qe(q) and write the full amplitude as 



A^\£) = drdyf T (r,y)M^ li (e,r,y) 



(D6) 



Results for M^ ull (£, r, y) can then be compared directly 
to the analogous results for — M EFT (r, y)S'(£) obtained 
with various M EFT results from App. [C] Here M FuI1 still 



depends on the thrust event shape because it contains 
6(£-Qe(pi)). 

The computation with one source and no additional 
gluons is very simple, and we find 



Full 



= -*(r,y) S'(£), 



(D7) 



where $(r, y) is given in Eq. (|C2|) . Next consider the 
computation for thrust at one loop. One performs the 
master-integral decomposition in the same way as for the 
anomalous dimension computation. The sum of all full 
theory virtual diagrams is UV finite and /^-independent 
(once one performs the usual QCD renormalization) . In 
fact the sum of diagrams involving a virtual gluon are 
identical in the full and EFT computations, and hence 
these contributions cancel when subtracting to determine 
the matching. 

Thus we only need master integrals involving full the- 
ory diagrams with real radiation to complete the calcu- 
lation. For these contributions only two new master in- 
tegrals are required to compute the /x-dependent pieces. 
For the thrust measurement on the real radiation gluon 
we will use the short-hand notation 



X T (£,p±) = - 5(£-p+)6(p--p + ) (D8) 
at L 



+ 6{£~p-)0(p+-p-) 



and to expand the master integrals with q <C £ we will 
use the identity 



S(x) ln(<J) 



ln(x) 



(D9) 



The first master integral shows up in the independent 
emission diagrams of Fig.|6]and when expanded for q <C £ 
gives 



~2e / d rf -y 



1 



1 



2\p\(2n) d - 1 p+ + Ap- +B 



(D10) 



-1 d 1 



(4tt) 2 d£ 



2 (p 
/' 



In 



AB 



7 ln l7 



S(£) 



2tt2 
3 



The second master integral appears in the double cut 
graphs involving the triple gluon vertex in Fig. and 
when expanded for q <C £ reads 



d d-i 



P 



1 



1 



2|p|(27r) rf - 1 p+ + A 2p ■ q + m 2 



M T (£,p ± ) 



-1 1_ d_ 

(4tt) 2 q+ dl' 



In 



7+2 



24 



?-(?) 



(Dll) 

are \i- independent 



where the omitted terms in the 4 
and UV-finite. 

Adding up all the /x-dependent full theory (real radia- 
tion) diagrams we find 



6'(£) lnf- 
V fx 



\- 


1 


M" 




\di 




.1. 


+ 



(D12) 



Note that \x d/d/i (r, y) = since the two contribu- 

tions cancel each other. This is consistent with the fact 



that there was no UV or IR divergence regulated by e. 
We still refer to them as /^-dependent terms since they 
have different functional dependence on I, In Eq. (|D12j) fx 
is simply a place holder scale for splitting the result into 
+-function and 5-function terms. For the corresponding 
EFT result, using Eq. (|UTT|) to get the sum of rcnormal- 
ized diagrams, we have 



'Hoop 



-S'(0(r, yfe^Ml - r 2 ) In ( -) 

7T V U J 



M ■ 

+ •••. (D13) 

When we subtract Eq. (|D13]> from Eq. ffD12]) we are left 
with only the /i-dependent +-function term. Identify- 
ing J dr dy f T (r,y)(fr(r,y) = this reproduces the /i- 
dependent term in the matching result given in Eq. (|7ip . 
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